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Abstract 

In this paper we study metastability in large volumes at low temperatures. We consider 
both Ising spins subject to Glauber spin-flip dynamics and lattice gas particles subject 
to Kawasaki hopping dynamics. Let j3 denote the inverse temperature and let C 1? 
be a square box with periodic boundary conditions such that lim^—^oo |A^| = oo. We run 
the dynamics on starting from a random initial configuration where all the droplets 
(= clusters of plus-spins, respectively clusters of particles) are small. For large /3, and for 
interaction parameters that correspond to the metastable regime, we investigate how the 
transition from the metastable state (with only small droplets) to the stable state (with 
one or more large droplets) takes place under the dynamics. This transition is triggered 
by the appearance of a single critical droplet somewhere in A^. Using potential-theoretic 
methods, we compute the average nucleation time {— the first time a critical droplet 
appears and starts growing) up to a multiplicative factor that tends to one as (5 — > oo. It 
turns out that this time grows as Ke T P for Glauber dynamics and K(3e T ^ for 
Kawasaki dynamics, where T is the local canonical, respectively, grand-canonical energy 
to create a critical droplet and if is a constant reflecting the geometry of the critical 
droplet, provided these times tend to infinity (which puts a growth restriction on lA^Q. 
The fact that the average nucleation time is inversely proportional to lA^I is referred to as 
homogeneous nucleation, because it says that the critical droplet for the transition appears 
essentially independently in small boxes that partition A^. 
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1 Introduction and main results 



1.1 Background 

In a recent series of papers, Gaudilliere, den Hollander, Nardi, Olivieri, and Scoppola [J_2j 
H3J E] study a system of lattice gas particles subject to Kawasaki hopping dynamics in a 
large box at low temperature and low density. Using the so-called path-wise approach to 
metastability (see Olivieri and Vares [23]), they show that the transition time between the 
metastable state (= the gas phase with only small droplets) and the stable state (= the liquid 
phase with one or more large droplets) is inversely proportional to the volume of the large 
box, provided the latter does not grow too fast with the inverse temperature. This type 
of behavior is called homogeneous nucleation, because it corresponds to the situation where 
the critical droplet triggering the nucleation appears essentially independently in small boxes 
that partition the large box. The nucleation time (= the first time a critical droplet appears 
and starts growing) is computed up to a multiplicative error that is small on the scale of 
the exponential of the inverse temperature. The techniques developed in [12} [T3l PH] center 
around the idea of approximating the low temperature and low density Kawasaki lattice gas 
by an ideal gas without interaction and showing that this ideal gas stays close to equilibrium 
while exchanging particles with droplets that are growing and shrinking. In this way, the 
large system is shown to behave essentially like the union of many small independent systems, 
leading to homogeneous nucleation. The proofs are long and complicated, but they provide 
considerable detail about the typical trajectory of the system prior to and shortly after the 
onset of nucleation. 

In the present paper we consider the same problem, both for Ising spins subject to Glauber 
spin-flip dynamics and for lattice gas particles subject to Kawasaki hopping dynamics. Using 
the potential-theoretic approach to metastability (see Bovier [5]), we improve part of the results 
in [1'2\ [TH1 [Ti] , namely, we compute the average nucleation time up to a multiplicative error 
that tends to one as the temperature tends to zero, thereby providing a very sharp estimate 
of the time at which the gas starts to condensate. 

We have no results about the typical time it takes for the system to grow a large droplet 
after the onset of nucleation. This is a hard problem that will be addressed in future work. 
All that we can prove is that the dynamics has a negligible probability to shrink down a su- 
percritical droplet once it has managed to create one. At least this shows that the appearance 
of a single critical droplet indeed represents the threshold for nucleation, as was shown in 
[12\ [P3"l E] • A further restriction is that we need to draw the initial configuration according 
to a class of initial distributions on the set of subcritical configurations, called the last-exit 
biased distributions, since these are particularly suitable for the use of potential theory. It 
remains a challenge to investigate to what extent this restriction can be relaxed. This problem 
is addressed with some success in [121 031 E] > an d will also be tackled in future work. 

Our results are an extension to large volumes of the results for small volumes obtained in 
Bovier and Manzo [8], respectively, Bovier, den Hollander, and Nardi [7j. In large volumes, 
even at low temperatures entropy is competing with energy, because the metastable state and 
the states that evolve from it under the dynamics have a highly non-trivial structure. Our 
main goal in the present paper is to extend the potential-theoretic approach to metastability 
in order to be able to deal with large volumes. This is part of a broader programme where 
the objective is to adapt the potential-theoretic approach to situations where entropy cannot 
be neglected. In the same direction, Bianchi, Bovier, and Ioffe [3] study the dynamics of the 
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random field Curie- Weiss model on a finite box at a fixed positive temperature. 

As we will see, the basic difficulty in estimating the nucleation time is to obtain sharp 
upper and lower bounds on capacities. Upper bounds follow from the Dirichlet variational 
principle, which represents a capacity as an infimum over a class of test functions. In [3] a 
new technique is developed, based on a variational principle due to Berman and Konsowa [2], 
which represent a capacity as a supremum over a class of unit flows. This technique allows 
for getting lower bounds and it will be exploited here too. 



1.2 Ising spins subject to Glauber dynamics 

We will study models in finite boxes, A^, in the limit as both the inverse temperature, (3, 
and the volume of the box, |A^|, tend to infinity. Specifically, we let C 1? be a square 
box with odd side length, centered at the origin with periodic boundary conditions. A spin 
configuration is denoted by a = {o~{x) : x G A^}, with a(x) representing the spin at site x, and 
is an element of Xa = {—1, +1} P. It will frequently be convenient to identify a configuration 
a with its support, defined as supp[<r] = {x G Ag: a(x) = +1}. 

The interaction is defined by the the usual Ising Hamiltonian 

Hp(°) = -{ E '(sMjO - \ E "W, vZXp, (1.1) 

x^y 

where J > is the pair potential, h > is the magnetic field, and x ~ y means that x and y 
are nearest neighbors. The Gibbs measure associated with Hg is 

^{a) = ^e~^\ aeXp, (1.2) 

where Zg is the normalizing partition function. 

The dynamics of the model will the a continuous-time Markov chain, (o~(t))t>o, with state 
space Xg whose transition rates are given by 

[ e -/3[H/3^')-H (a)] + for a' = a x for some x £ Ag, M 
C ^ a) = \0, otherwise, (L3) 

where a x is the configuration obtained from a by flipping the spin at site x. We refer to this 
Markov process as Glauber dynamics. It is ergodic and reversible with respect to its unique 
invariant measure, fig, i.e., 

lJ-p{a)cp(o-,a') = ng{a')cg{a \a), \/<T,a' € Xg. (1.4) 



Glauber dynamics exhibits metastable behavior in the regime 

< h < 2J, (3 oo. (1.5) 

To understand this, let us briefly recall what happens in a finite /3-independent box A C Z 2 , 
Let Ba and EBa denote the configurations where all spins in A are —1, respectively, +1. As 
was shown by Neves and Schonmann [22], for Glauber dynamics restricted to A with periodic 
boundary conditions and subject to (jl.5p . the critical droplets for the crossover from Ba to Ba 
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Figure 1: A critical droplet for Glauber dynamics on A. 
the non-shaded area the (— l)-spins (see (|1.6jl ). 



The shaded area represents the (+l)-spins, 



are the set of all those configurations where the (+l)-spins form an £ c x [l c — 1) quasi-square 
(in either of both orientations) with a protuberance attached to one of its longest sides, where 



2J 
~~h 



(1.6) 



(see Figs.[Hand[2l for non-degeneracy reasons it is assumed that 2J/h £ N). The quasi-squares 
without the protuberance are called proto- critical droplets. 

Let us now return to our setting with finite /3-dependent volumes Ap C Z 2 . We will start 
our dynamics on from initial configurations in which all droplets are "sufficiently small" . To 
make this notion precise, let Cb(o~), a G Xp, be the configuration that is obtained from a by a 
"bootstrap percolation map", i.e., by circumscribing all the droplets in a with rectangles, and 
continuing to doing so in an iterative manner until a union of disjoint rectangles is obtained 
(see Kotecky and Olivieri [H]). We call Cb(o~) subcritical if all its rectangles fit inside a 
proto-critical droplet and are at distance > 2 from each other (i.e., are non-interacting). 



Definition 1.1 (a) S = {a G Xp: Cb(o~) is subcritical}. 

(b) V = {a G S : cp(a,a') > for some a' G S c }. 

(c) C = {a' G S c : cp(a,a') > for some a G S}. 



We refer to S, V and C as the set of subcritical, proto-critical, respectively, critical configu- 
rations. Note that, for ever a G Xp, each step in the bootstrap percolation map a — > Cb(<j) 
deceases the energy, and therefore the Glauber dynamics moves from a to Cb{o~) in a time of 
order one. This is why Cb(o~) rather than a appears in the definition of S. 

For £i,£2 G N, let Ri 1: £ 2 (x) C A^ be the l\ x £2 rectangle whose lower-left corner is x. We 
always take l\ < £2 and allow for both orientations of the rectangle. For L = 1, . . . , 2l c — 3, let 
Ql{ x ) denote the L-th element in the canonical sequence of growing squares and quasi-squares 

Rl,2(x), R 2 ,2(x), R2,3(x), #3,3 (a?), ■ • • > #4-1,4-1 ( X )> #4-1,40*0- ( L7 ) 

In what follows we will choose to start the dynamics in a way that is suitable for the use of 
potential theory, as follows. First, we take the initial law to be concentrated on sets Sl C S 
defined by 

Sl = {o~ G S: each rectangle in Cb(ct) fits inside Ql(x) for some x G A^} , (1-8) 
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where L is any integer satisfying 



L* < L < 2£ c - 3 with L* = min <l<L<2£ c — 3: lim Lj = 1 I . (1.9) 



In words, Sl is the subset of those subcritical configurations whose droplets fit inside a square 
or quasi-square labeled L, with L chosen large enough so that Sl is typical within S under 
the Gibbs measure \ip as (3 — > 00 (our results will not depend on the choice of L subject to 
these restrictions). Second, we take the initial law to be biased according to the last exit of 
Sl for the transition from Sl to a target set in S c . (Different choices will be made for the 
target set, and the precise definition of the biased law will be given in Section F2.21 ) This is a 
highly specific choice, but clearly one of physical interest. 

Remarks: (1) Note that S2e c -3 = S, which implies that the range of L-values in (|1.9p 
is non-empty. The value of L* depends on how fast Ap grows with (3. In Appendix IC.ll 
we will show that, for every 1 < L < 2£ c — 4, lim^oo ^lq{Sl) / ^g{S) = 1 if and only if 
lim^oo \kg\e~ l3VL + 1 = with T^+i the energy needed to create a droplet Ql+i(0) at the 
origin. Thus, if (A^ = e e ^, then L* = L*(9) = (24 - 3) A min{L G N: T L+ i > 0}, which 
increases stepwise from 1 to 2£ c — 3 as 6 increases from to T defined in (jl.lOj) . 
(2) If we draw the initial configuration <to from some subset of S that has a strong recurrence 
property under the dynamics, then the choice of initial distribution on this subset should not 
matter. This issue will be addressed in future work. 



Figure 2: A nucleation path from El a to EBa for Glauber dynamics. T in (ll.lOj) is the minimal energy 
barrier the path has to overcome under the local variant of the Hamiltonian in 

To state our main theorem for Glauber dynamics, we need some further notation. The 
key quantity for the nucleation process is 



which is the energy needed to create a critical droplet of (+l)-spins at a given location in 
a sea of (— l)-spins (see Figs. [T]and[5]). For a G Xg, let P CT denote the law of the dynamics 
starting from a and, for v a probability distribution on X, put 




r = j[u c ] - M4(4 - 1) + 1] 



(1.10) 




(1.11) 



For a non-empty set A C Xg, let 



t a = inf{i > 0: a t G A, a t - £ A} 



(1.12) 
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denote the first time the dynamics enters A. For non-empty and disjoint sets A,B<Z Xg, 
let Vj^ denote the last-exit biased distribution on A for the crossover to B defined in (|2.9p in 
Section O Put 

N X =U C , N 2 = |(24-1). (1.13) 
For M £ N with M > £ c , define 

T>m = {o- G Xp: 3x£ A/3 such that supp[C B (cr)] D Rm,m{x)}, (1-14) 

i.e., the set of configurations containing a supercritical droplet of size M. For our results 
below to be valid we need to assume that 

lim \K P \ = oo, lim |Afl| e _/3r = 0. (1.15) 

/3— »oo /3^oo 

Theorem 1.2 7n the regime (|1.5p . subject to (|1.9p and (jl. 15[) . i/ie following hold: 
(a) 

lim lA^e-^E- ( r<sc ) = -L. (1.16) 



lim |A^| e-^E S c XC (r^) = -L (1.17) 

p^oo "s L iV2 



lim |A^| e-^E OM (rp M ) = V£ c < M < 24 - 1. (1.18) 

/3^oo Ks L iV2 



The proof of Theorem 11.21 will be given in Section [3l Part (a) says that the average time 
to create a critical droplet is [1 + o(l)]e^ r /A r i|A j a|. Parts (b) and (c) say that the average 
time to go beyond this critical droplet and to grow a droplet that is twice as large is [1 + 
o(l)]e^ r /A r 2|A ( g|. The factor N± counts the number of shapes of the critical droplet, while 
lA^I counts the number of locations. The average times to create a critical, respectively, a 
supercritical droplet differ by a factor N2/N1 < 1. This is because once the dynamics is "on 
top of the hill" C it has a positive probability to "fall back" to S. On average the dynamics 
makes N1/N2 > 1 attempts to reach the top C before it finally "falls over" to S C \C. After 
that, it rapidly grows a large droplet (see Fig. [2]). 

Remarks: (1) The second condition in (|1.15p will not actually be used in the proof of 
Theorem 11.2( a). If this condition fails, then there is a positive probability to see a proto- 
critical droplet in A^ under the starting measure v§ , and nucleation sets in immediately. 
Theorem 11.2( a) continues to be true, but it no longer describes metastable behavior. 

(2) In Appendix [D] we will show that the average probability under the Gibbs measure \ip 
of destroying a supercritical droplet and returning to a configuration in Sl is exponentially 
small in (3. Hence, the crossover from Sl to S C \C represents the true threshold for nucleation, 
and Theorem ll.2f b) represents the true nucleation time. 

(3) We expect Theorem 11.2( c) to hold for values of M that grow with j3 as M = e ^' . 
As we will see in Section 13.31 the necessary capacity estimates carry over, but the necessary 
equilibrium potential estimates are not yet available. This problem will be addressed in future 
work. 

(4) Theorem 11.21 should be compared with the results in Bovier and Manzo [8] for the case 
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of a finite /3-independent box A (large enough to accommodate a critical droplet). In that 
case, if the dynamics starts from Ba, then the average time it needs to hit C\ (= the set of 
configurations in A with a critical droplet), respectively, Ba equals 

Ke^[l + o(l)), with K = K(A,l c ) = 1 JL for N = N U N 2 . (1.19) 

(4) Note that in Theorem 11.21 we compute the first time when a critical droplet appears 
anywhere (!) in the box A«. It is a different issue to compute the first time when the plus- 
phase appears near the origin. This time, which depends on how a supercritical droplet grows 
and eventually invades the origin, was studied by Dehghanpour and Schonmann [101 lllj . 
Shlosman and Schonmann |24| and, more recently, by Cerf and Manzo [9]. 



1.3 Lattice gas subject to Kawasaki dynamics 

We next consider the lattice gas subject to Kawasaki dynamics and state a similar result for 
homogeneous nucleation. Some aspects are similar as for Glauber dynamics, but there are 
notable differences. 

A lattice gas configuration is denoted by a = {a{x) : x £ Xp}, with a(x) representing the 
number of particles at site x, and is an element of Xp = {0, l} 7 ^- 3 . The Hamiltonian is given 
by 

Hp{<r) = -U Y, °(, x Xv)> <reXp, (1-20) 

where —U < is the binding energy and x ~ y means that x and y are neighboring sites. 
Thus, we are working in the canonical ensemble, i.e., there is no term analogous to the second 
term in (|1.1|) . The number of particles in A^ is 

np=\pp\A p \l (1.21) 

where pp is the particle density, which is chosen to be 

pp = e _/3A , A > 0. (1.22) 

Put 

Xf^ = {aeXp: |suppH| = n^}, (1.23) 
where supp[<r] = {x £ Ap : a(x) = 1}. 

Remark: If we were to work in the grand- canonical ensemble, then we would have to consider 
the Hamiltonian 



with A > an activity parameter taking over the role of h in (|1.1|) . The second term would 
mimic the presence of an infinite gas reservoir with density pp outside A^. Such a Hamiltonian 
was used in earlier work on Kawasaki dynamics, when a finite /3-independent box with open 
boundaries was considered (see e.g. den Hollander, Olivieri, and Scoppola [IB], den Hollander, 
Nardi, Olivieri, and Scoppola [T7], and Bovier, den Hollander, and Nardi [7]). 



7 



The dynamics of the model will be the continuous-time Markov chain, (<7t)t>o, with state 
space X^ 1 ^ whose transition rates are 

, i \ f e -i3[^K)-^Wl+ for a' = a x ' y for some x, y £ K with x ~ y, ot ,, 
C ^ a) = \0, otherwise, (L25) 

where a x ' y is the configuration obtained from a by interchanging the values at sites x and 
y. We refer to this Markov process as Kawasaki dynamics. It is ergodic and reversible with 
respect to the canonical Gibbs measure 

p p (<r) = -^e-f> B »M aeX^\ (1.26) 

where zj^ 13 ^ is the normalizing partition function. Note that the dynamics preserves particles, 
i.e., it is conservative. 




A 

Figure 3: A critical droplet for Kawasaki dynamics on A (— a proto-critical droplet plus a free particle). 
The shaded area represents the particles, the non-shaded area the vacancies (see (|1.28[0 . Note that 
the shape of the proto-critical droplet for Kawasaki dynamics is the same as that of the critical droplet 
for Glauber dynamics. The proto-critical droplet for Kawasaki dynamics becomes critical when a free 
particle is added. 

Kawasaki dynamics exhibits metastable behavior in the regime 

U < A < 2U, P^oo. (1.27) 

This is again inferred from the behavior of the model in a finite /3-independent box A C Z 2 . Let 
□a and M\ denote the configurations where all the sites in A are vacant, respectively, occupied. 
For Kawasaki dynamics on A with an open boundary, where particles are annihilated at rate 
1 and created at rate e _A ^, it was shown in den Hollander, Olivieri, and Scoppola [18] and in 
Bovier, den Hollander, and Nardi [7j that, subject to (jl.27|) and for the Hamiltonian in (jl.24p . 
the critical droplets for the crossover from H\\ to M\ are the set of all those configurations 
where the particles form 

(1) either an (£ c — 2) x (£ c — 2) square with four bars attached to the four sides with total 
length 3£ c — 3, 

(2) or an {l c — 1) x (£ c — 3) rectangle with four bars attached to the four sides with total 
length 34-2, 



S 



plus a free particle anywhere in the box, where 

U 



2U - A 



(1.28) 



(see Figs. [3] and HI for non-degeneracy reasons it is assumed that U/(2U — A) ^ N). 

Let us now return to our setting with finite /3-dependent volumes. We define a reference 
distance, Lp, as 

L 2 e (A- W = J_ e -^ (L29) 

with 

lim dp = 0, lim f38p = oo, (1.30) 

i.e., is marginally below the typical interparticle distance. We assume to be odd, and 
write BLg t Lp(x), x G A^, for the square box with side length Lp whose center is x. 

Definition 1.3 (a) S = {a G xf &) : |supp[<r] n B Lp>Lfl (x)\ < 4(4 - 1) + 1 V x G A/?}. 
(b) V = {a £ S : Cp{a,o J ) > /or some </ G 5 C }. 
fc,) C = {(j'6 5 c : cp(a,o J ) > /or some a G 5}. 

fdj C~ = {a G C: 3x G A^ suc/i that Bl^^lAx) contains a proto- critical droplet plus a free 
particle at distance Lp}. 

(e) C + = the set of configurations obtained from C~ by moving the free particle to a site at 
distance 2 from the proto- critical droplet. 



As before, we refer to S, V and C as the set of subcritical, proto- critical, respectively, critical 
configurations. Note that, for every a G S, the number of particles in a box of size Lp does 
not exceed the number of particles in a proto-critical droplet. These particles do not have 
to form a cluster or to be near to each other, because the Kawasaki dynamics brings them 
together in a time of order L\ = o(l/pp). 

The initial law will again be concentrated on sets Sl C S, this time defined by 

S L = {o- G Xf !i) : |supp[o-] n B LfhLp {x)\ < L Vx G A^}, (1.31) 
and L any integer satisfying 

L* < L < 4(4 - 1) + 1 with V = min ( 1 < L < 4(4 - 1) + 1 : lim ^f^) = ll . 

I P-+oo np{S) J 

(1.32) 

In words, Sl is the subset of those subcritical configurations for which no box of size Lp carries 
more than L particles, with L again chosen such that Sl is typical within S under the Gibbs 
measure up as (3 — > oo. 

Remark: Note that Sp c ^ c _i^ + i = S. As for Glauber, the value of L* depends on how 
fast Ap grows with f3. In Appendix IC2I we will show that, for every 1 < L < 4(4 — l)j 
lim^oo np(S L )/np(S) = 1 if and only if lfing-^ | A /3 |e~ /3(rL + 1_A) = with r L+ i the energy 
needed to create a droplet of L + 1 particles, closest in shape to a square or quasi-square, 
in BL p ,Lp(Q) under the grand-canonical Hamiltonian on this box. Thus, if [A^l = e 9 @, then 
L* = L*{9) = [4(4 — 1) + 1] A min{L G N: r^+i — A > 9}, which increases stepwise from 1 
to 4(4 — 1) + 1 as 9 increases from A to V defined in (jl.33p . 
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Figure 4: A nucleation path from 0\ to M\ for Kawasaki dynamics on A with open boundary, 
r in (|1.33[) is the minimal energy barrier the path has to overcome under the local variant of the 
grand-canonical Hamiltonian in (|1.24|) . 



Set 

r = -U[(£ c - l) 2 + 4(4 - 1) + 1] + A[4(4 - 1) + 2], (1.33) 

which is the energy of a critical droplet at a given location with respect to the grand- canonical 
Hamiltonian given by ([04]) (see Figs. EJandH]). Put N = \l 2 c {l 2 c - 1). For M € N with 
M > 4, define 

Dm = {cr £ X@: 3x G A/3 such that supp[(a)] D Rm,m{x)}, (1-34) 

i.e., the set of configurations containing a supercritical droplet of size M. For our results 
below to be valid we need to assume that 

lim | A/3 1 p/3 = do, lim (A^ e _/3r = 0. (1.35) 

13— >oo /3— >oo 

This first condition says that the number of particles tends to infinity, and ensures that the 
formation of a critical droplet somewhere does not globally deplete the surrounding gas. 

Theorem 1.4 In the regime (|1.27|) . subject to (ll.32|) and (II. 35ft . i/ie following hold: 

lApl JK e ~^ E 4^>c+ (^\c)uc+) = ^ (1-36) 

(b) 

lim |A^| e-^E v M (t Vm ) = 1 V4 < M < 24 - 1. (1-37) 

The proof of Theorem 1 1 .4^ which is the analog of Theorem 11.2} will be given in Section 21 
Part (a) says that the average time to create a critical droplet is [1 + o(l)](/3A/4-7r)e /3r iV|A / g|. 
The factor /3A/47T comes from the simple random walk that is performed by the free particle 
"from the gas to the proto-critical droplet" (i.e., the dynamics goes from C~ to C + ), which 
slows down the nucleation. The factor N counts the number of shapes of the proto-critical 



10 



droplet (see Bovier, den Hollander, and Nardi [7J). Part (b) says that, once the critical droplet 
is created, it rapidly grows to a droplet that has twice the size. 

Remarks: (1) As for Theorem ll.2l fc). we expect Theorem II. 4( b) to hold for values of M that 
grow with as M = . See Section O for more details. 

(2) In Appendix [D] we will show that the average probability under the Gibbs measure ug of 
destroying a supercritical droplet and returning to a configuration in Sl is exponentially small 
in f3. Hence, the crossover from Sl to S C \C U C + represents the true threshold for nucleation, 
and Theorem 11.4( a) represents the true nucleation time. 

(3) It was shown in Bovier, den Hollander, and Nardi [7j that the average crossover time in a 
finite box A equals 

KeF[l + o{l)), w ithir = ir(A,4)~^|^^pT, A^Z 2 . (1.38) 

This matches the |A^| -dependence in Theorem 11,41 with the logarithmic factor in (|1.38|) ac- 
counting for the extra factor /3A in Theorem 11.41 compared to Theorem 11.21 Note that this 
factor is particularly interesting, since it says that the effective box size responsible for the 
formation of a critical droplet is Lp. 



1.4 Outline 

The remainder of this paper is organized as follows. In Section [2] we present a brief sketch of 
the basic ingredients of the potential-theoretic approach to met ast ability. In particular, we 
exhibit a relation between average crossover times and capacities, and we state two variational 
representations for capacities, the first of which is suitable for deriving upper bounds and the 
second for deriving lower bounds. Section [3] contains the proof of our results for the case of 
Glauber dynamics. This will be technically relatively easy, and will give a first flavor of how our 
method works. In Section 0] we deal with Kawasaki dynamics. Here we will encounter several 
rather more difficult issues, all coming from the fact that Kawasaki dynamics is conservative. 
The first is to understand why the constant T, representing the local energetic cost to create 
a critical droplet, involves the grand-canonical Hamiltonian, even though we are working in 
the canonical ensemble. This mystery will, of course, be resolved by the observation that the 
formation of a critical droplet reduces the entropy of the system: the precise computation of 
this entropy loss yields T via equivalence of ensembles. The second problem is to control the 
probability of a particle moving from the gas to the proto-critical droplet at the last stage 
of the nucleation. This non-locality issue will be dealt with via upper and lower estimates. 
Appendices lAUDl collect some technical lemmas that are needed in Sections [3] [U 

The extension of our results to higher dimensions is limited only by the combinatorial 
problems involved in the computation of the number of critical droplets (which is hard in the 
case of Kawasaki dynamics) and of the probability for simple random walk to hit a critical 
droplet of a given shape when coming from far. We will not pursue this generalization here. 
The relevant results on a /3-independent box in Z 3 can be found in Ben Arous and Cerf pQ 
(Glauber) and den Hollander, Nardi, Olivieri, and Scoppola [T7] (Kawasaki). For recent 
overviews on droplet growth in metastability, we refer the reader to den Hollander [151 PTS] 
and Bovier [HE]. A general overview on metastability is given in the monograph by Olivieri 
and Vares [23]. 
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2 Basic ingredients of the potential-theoretic approach 



The proof of Theorems 11.21 and 11.41 uses the potential-theoretic approach to met ast ability 
developed in Bovier, Eckhoff, Gayrard and Klein [6] . This approach is based on the following 
three observations. First, most quantities of physical interest can be represented in term of 
Dirichlet problems associated with the generator of the dynamics. Second, the Green function 
of the dynamics can be expressed in terms of capacities and equilibrium potentials. Third, 
capacities satisfy variational principles that allow for obtaining upper and lower bounds in a 
flexible way. We will see that in the current setting the implementation of these observations 
provides very sharp results. 



2.1 Equilibrium potential and capacity 



The fundamental quantity in the theory is the equilibrium potential, h A ^ , associated with two 

non-empty 
is given by 



non-empty disjoint sets of configurations, A, B C X (= Xp or X^ 1 ^), which probabilistically 



P<x(ta<tb), for a € (A U Bf , 

h A ,B(a) = { 1, for o-EA, (2.1) 

0, for a £ B, 

where 

t a = mf{t > 0: a t eA,a t - <£ A}, (2.2) 

(o~t)t>o is the continuous-time Markov chain with state space X, and F a is its law starting 
from a. This function is harmonic and is the unique solution of the Dirichlet problem 

(Lh A , B )(a) = 0, a£(AuB) c , 

h A ,B{°) = ffU (2.3) 
hAjsi?) = °> cr e B, 

where the generator is the matrix with entries 

L(cr, a') = cp(a, a') - 5 a ^ cp(cr), a, a' £ X, (2.4) 
with cp(a) the total rate at which the dynamics leaves a, 

c p(°)= Yl c /3(^')> a ^ x - ( 2 - 5 ) 

a'eX\{a} 

A related quantity is the equilibrium measure on ^4, which is defined as 

e A ,B^) = -{Lh AB ){a), a e A. (2.6) 
The equilibrium measure also has a probabilistic meaning, namely, 

K(tb<t a ) = ^^-, aeA. (2.7) 

The key object we will work with is the capacity, which is defined as 

CAP(A, B) = J2 ^(ff)^s(^). (2-8) 
<re A 
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2.2 Relation between crossover time and capacity 



The first important ingredient of the potential-theoretic approach to metastability is a formula 
for the average crossover time from A to B. To state this formula, we define the probability 
measure Vj^ on A we already referred to in Section [TJ namely, 



CAP(AB) . for a G A 
0, for cr <E .A c . 

The following proposition is proved e.g. in Bovier [5]. 
Proposition 2.1 For any two non-empty disjoint sets A, B C X, 



2.i 



y%y) ^{t b ) = J, B) Yl ^(*) h A,B(°)- (2-io) 

Remarks: (1) Due to (l2.7H2.8p . the probability measure z^(c) can be written as 

^) = ^^^B<r A ), a e A, (2.11) 

and thus has the flavor of a last-exit biased distribution. Proposition 12.11 explains why our 
main results on average crossover times stated in Theorem 1 1 . 2 1 and 1 1 . 41 are formulated for this 
initial distribution. Note that 



<T&B C 

We will see that in our setting /j,^(B c \A) = o([j,p(A)) as (3 — > oo, so that the sum in the 
right-hand side of (|2,10p is ~ fig (A) and the computation of the crossover time reduces to the 
estimation of CAP(A,B). 

(2) For a fixed target set B, the choice of the starting set A is free. It is tempting to choose 
A = {cr} for some cr G X. This was done for the case of a finite /3-independent box A. 
However, in our case (and more generally in cases where the state space is large) such a choice 
would give intractable numerators and denominators in the right-hand side of (|2.1U|) . As a 
rule, to make use of the identity in f)2. 10p . A must be so large that the harmonic function h^B 
"does not change abruptly near the boundary of A" for the target set B under consideration. 

As noted above, average crossover times are essentially governed by capacities. The use- 
fulness of this observation comes from the computability of capacities, as will be explained 
next. 



2.3 The Dirichlet principle: A variational principle for upper bounds 

The capacity is a boundary quantity, because e^,B > only on the boundary of A. The 
analog of Green's identity relates it to a bulk quantity. Indeed, in terms of the Dirichlet form 
defined by 

£{h) = \ J2 ^)cB{a,a')[h{a)-h{a')}\ h:X^[0,l], (2.13) 

a,cr'eX 
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it follows, via (JET} and Pig} , that 

CAP(AS) = 5(/i AB ). (2.14) 
Elementary variational calculus shows that the capacity satisfies the Dirichlet principle: 



Proposition 2.2 For any too non-empty disjoint sets A, B C X, 

CAP{A,B) = min (2.15) 

fc: #->[0,l] 

The importance of the Dirichlet principle is that it yields computable upper bounds for capaci- 
ties by suitable choices of the test function h. In metastable systems, with the proper physical 
insight it is often possible to guess a reasonable test function. In our setting this will be seen 
to be relatively easy. 



2.4 The Berman-Konsowa principle: A variational principle for lower 
bounds 

We will describe a little-known variational principle for capacities that is originally due to 
Berman and Konsowa [2]. Our presentation will follow the argument given in Bianchi, Bovier, 
and Ioffe PJ. 

In the following it will be convenient to think of X as the vertex-set of a graph (X,£) 
whose edge-set £ consists of all pairs (a, a'), a, a' E X, for which cp{p,o') > 0. 

Definition 2.3 Given two non-empty disjoint sets A, B C X, a loop-free non-negative unit 
flow, f , from A to B is a function f : £ — * [0, oo) such that: 
(a;(/(e)>0=>/(-e) = 0)Ve€£. 

(b) f satisfies Kirchoff's law: 

f(o-,o-')= /(*'», V<TEX\(AUB). (2.16) 

(c) f is normalized: 

EE/(^) = 1= EE/(^ ff )- ( 2 - 17 ) 

(d) Any path from A to B along edges e such that /(e) > is self-avoiding. 

The space of all loop-free non-negative unit flows from A to B is denoted by U^g. 

A natural flow is the harmonic flow, which is constructed from the equilibrium potential 
h AtB as 

fAJBfa o-') = CAp ^ A ff) Pfi(<r)cp(<r, o-') [hAjsia) - h AB (a')] + , a, a' G X. (2.18) 

It is easy to verify that /^g satisfies (a-d). Indeed, (a) is obvious, (b) uses the harmonicity 
of h^B, (c) follows from (|2.6j) and (|2.8p . while (d) comes from the fact that the harmonic flow 
only moves in directions where Hj^b decreases. 
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A loop-free non-negative unit flow / is naturally associated with a probability measure 
P* on self-avoiding paths, 7. To see this, define F(a) = Yla'ex f( a i a ')> a ^ X\B. Then 

is the Markov chain (o" n )neN with initial distribution P^(o"o) = F(gq)\a{gq), transition 
probabilities 



q f (a, a') 



F(a) ' 



a E X\B, 



(2.19) 



such that the chain is stopped upon arrival in B. In terms of this probability measure, we 
have the following proposition (see [3] for a proof). 

Proposition 2.4 Let A,B <Z X be two non-empty disjoint sets. Then, with the notation 
introduced above, 



CAP(A, B) = sup W 



_e€7 



f(ei,e r ) 



V0(ei)cp(ei,e r 



(2.20) 



where e = (ej, e r ) and the expectation is with respect to 7. Moreover, the supremum is realized 
for the harmonic flow Ja,b- 



The nice feature of this variational principle is that any flow gives a computable lower 
bound. In this sense (I2.15D and (|2.20p complement each other. Moreover, since the harmonic 
flow is optimal, a good approximation of the harmonic function h^s by a test function h leads 
to a good approximation of the harmonic flow f^s by a test flow / after putting h instead 
of n A,B i n (|2.18p . Again, in metastable systems, with the proper physical insight it is often 
possible to guess a reasonable flow. We will see in Sections [3HH how this is put to work in our 
setting. 

3 Proof of Theorem 11.2 
3.1 Proof of Theorem OK a) 

To estimate the average crossover time from Sl C S to S c , we will use Proposition 12.11 With 
A = S L and B = S c , (|2~10l) reads 

J~ ug(a)E a (r S c) = 1 — ^2^{a)hs L ^(a). (3.1) 

The left-hand side is the quantity of interest in (j!.16[) . In Sections I3.1.1I[3~1.2I we estimate 
So-g5 ^/3(°")^'5l,5 c ( (J ) an d CAP(Sl,S c ). The estimates will show that 

r.h.s. TO = — e^ r [1 + 0(1)], (3.2) 
J\l\A/3\ 

3.1.1 Estimate of Y^aeS ^0{ a ) h s L ,s^) 

Lemma 3.1 Y,aes W( a ) h S L ,S°((r) = Pp{S)[l + o(l)\ as (3 -> 00. 
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Proof. Write, using (|2.ip . 

^2 l ^{a)hs L ,s-{o-)=^2 l Hp{a-)hs L ,s-{o-)+ Y fj,p(a)hs L ,s-((r) 

aes aes L &es\s L ^ 

= Hp{S L ) + Y Vp(?Wv{ts l < r S c). 
<?es\s L 

The last sum is bounded above by fj,p(S\Sz). But fj,p(S\Si) = o(fip(S)) as (3 — > oo by our 
choice of L in (11.91) . I 



3.1.2 Estimate of CAP(S L ,S C ) 

Lemma 3.2 CAP(S L ,S C ) = N x |A /3 | e - /3r /x /3 (5)[l + o(l)] as /? -» oo toitfi JVi = 44. 
Proof. The proof proceeds via upper and lower bounds. 

Upper bound: We use the Dirichlet principle and a test function that is equal to 1 on S to 
get the upper bound 

cap(s l ,s c ) < cap(s,s c ) = Y M^M^V) = Y M*) A/i/sO/)] ^ M c )< 

o-es,<T'GS c o-es : o-'es c 

o | g(ir,CT')>0 c^(ct,ct')>0 

(3.4) 

where the second equality uses (I1.4D in combination with the fact that cp(a, a') V cp(a', a) = 1 
by (|1.3p . Thus, it suffices to show that 

A*/3(C) < A f i|A /3 |e- /3r [l + o(l)] as (3 oo. (3.5) 

For every a £V there are one or more rectangles Ri c -i t e c (x), x = x(a) £ Xp, that are filled by 
(+l)-spins in Cb{o~). If o' E C is such that a' = a y for some y £ Ap, then a' has a (+l)-spin 
at y situated on the boundary of one of these rectangles. Let 

S(x) = {a £S: supp[a] C R ic _ 1A (x)}, 

S(x) = {a£S: suppH C [R ic+1/c+2 (x - (l,l))] c }. 




Figure 5: Rt e —i t i e (x) (shaded box) and [Ri c+ u c+ 2(x — (1, l))] c (complement of dotted box). 

For every a £ V, we have a = a V a for some a £ <S(x) and a £ S(x), uniquely decomposing 
the configuration into two non-interacting parts inside Ri c -i ) i c (x) and [Rt c+ i^ c+ 2 (x — (1, l))] c 
(see Fig. [S]) . We have 

Hp(a) - Hp(B) = [Hp(a) - Hp(B)] + [Hp(a) - Hp(B)]. (3.7) 
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Moreover, for any y £ supp[Cs (a)], we have 



Hp{a y )>Hp(a) + 2J-h. (3.8) 

Hence 

mo = j- E E ^"» ( -> 

a*eC 

* leAf, d- Sl s(x) «es(x) (3.9) 

< [i + (i)] 1 ^lA^e"^ V e-mm 

p <reS(o) 
= [l + o(l)]JV 1 |A /9 |e-^ r ^(0)), 

where the first inequality uses ()3.7H3.8p . with Aq = 2 x 2£ c = 4£ c counting the number of 
critical droplets that can arise from a proto-critical droplet via a spin flip (see Fig. [1]), and 
the second inequality uses that 

a e 5(0), ffVffGP^ H p (a) > H p (Re c -i /c (0)) = T — (2 J — h) + Hp(B) (3.10) 

with equality in the right-hand side if and only if supp[<r] = Ri c -i £ c (0). Combining (|3.4p and 
(|3.9j) with the inclusion 5(0) C 5, we get the upper bound in (13. 5p . 

Lower bound : We exploit Proposition 12.41 by making a judicious choice for the flow /. In fact, 
in the Glauber case this choice will be simple: with each configuration a £ Sl we associate 
a configuration in C C 5 C with a unique critical droplet and a flow that, from each such 
configuration, follows a unique deterministic path along which this droplet is broken down 
in the canonical order (see Fig. |SJ) until the set Sl is reached, i.e., a square or quasi-square 
droplet with label L is left over (recall (jl.7H1.8p ). 



Ql 



05 <*K 

Figure 6: Canonical order to break down a critical droplet. 
Let f{j3) be such that 

hm/(/3) = oo, lim - log /(/?) = 0, lim |A^| //(/?) = oo, (3.11) 

p^oo p— >oo p p— >oc 

and define 

W = {a e 5: |suppH| < |A/3|//(/3)}. (3.12) 

Let Cl C C C 5 C be the set of configurations obtained by picking any a £ Sl n W and 
adding somewhere in A/j a critical droplet at distance > 2 from supp[er]. Note that the 
density restriction imposed on W guarantees that adding such a droplet is possible almost 
everywhere in Ap for (5 large enough. Denoting by Pr y \{x) the critical droplet obtained by 
adding a protuberance at y along the longest side of the rectangle g c (x), we may write 

C L = {a U P (y) (x) : a G 5 n W, x,y E Ap, (x,y)J_<r}, (3.13) 
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where (x,y)-La stands for the restriction that the critical droplet Pr y \(x) is not interacting 
with supp[<r], which implies that Hp (a U Pi y \{x)) = Hp{a) + T (see Figs. [7]and[H|). 



Figure 7: The critical droplet P( y ){x 



□ 



Figure 8: Going from Sl to Cl by adding a critical droplet P( y )(x) somewhere in A^. 

Now, for each a G Cl, we let 7<j = (7^(0), 7 CT (1), . . . , ~f a (K)) be the canonical path from 
a = 7cr(0) to Sl along which the critical droplet is broken down, where K = v(2£ c — 3) — v{L) 
with 

v(L) = \Q L (0)\ (3.14) 
(recall <|1.7j) ). We will choose our flow such that 

fW,G") 

vq(g), if g' = g, g" = 7o- (1) for some g G Cl, 

SseC L fila(k - l),~/ a (k)), if g' = 7 <T (fc), g" = 7a- (k + l)for some k > 1, g G Cl, 
0, otherwise. 

(3.15) 

Here, fo is some initial distribution on Cl that will turn out to be arbitrary as long as its 
support is all of Cl- 

We see from (I3.15D that the flow increases whenever paths merge. In our case this happens 
only after the first step, when the protuberance at y is removed. Therefore we get the explicit 
form 

fo(cr), if g 1 = g, a" = 7ct(1) for some g G Cl, 
/(g', g") = { Cv (g), ifcr' = 7 <T (Jfe), g" = 7 ff (fc + 1) for some k > 1, g € C L , (3.16) 
0, otherwise, 
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where C = 2£ c is the number of possible positions of the protuberance on the proto-critical 
droplet (see Fig. [H]) . Using Proposition 12.41 we therefore have 



cap(S l ,S c ) = cap(S c ,S l ) > cap(C l ,S l ) 

[K-\ 



°ec L 



E 



/(7<r(k),7»(* + l)) 



-1 



E 



~ M/9(7«r(fc))Ci9(7«T(fe)i7ff(fc+l)) 
1 K ~ 1 



(3.17) 

Thus, all we have to do is to control the sum between square brackets. 

Because c i g(7 CT (0),7 o -(l)) = 1 (removing the protuberance lowers the energy), the term 
with k = equals 1/ np(cr). To show that the terms with k > 1 are of higher order, we argue 
as follows. Abbreviate S = /i(^ c — 2). For every k > 1 and <r(0) G Cl, we have (see Fig. [9] and 
recall IjOOl) ) 



^(7,(^)^(7,(^,7.(^+1)) = ^- e -m^))vH fjMk+m > Mao)e mj-h-s] 

where c5 = 2J-/i-~=2J- h(£ c - 1) > (recall fity ). Therefore 



E 



C 



< 



1 



fc J^ fJ,^ a (k))c l3 (^ (T (k), J<r(k + 1)) ^((T) 

and so from (13.171) we get 



CKe 



-8(3 



(3.18) 
(3.19) 



CAP(5 L ,5 C )> Y, J 



<rec L 



+ CKe-l 35 1 + CKe~P & 



[l + o(l)]^(C L ) 



(3.20) 




Figure 9: Visualization of (|3.18D . 
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The last step is to estimate, with the help of (|3.13j) . 

Zr — ' Zn — ' 

1 a£C L 1 crG5 L nW a.S/SA^ 

(x,y)±tr 



<~ pv ^- y e-wfa) y 1 ( 3 - 21 ) 



(x,y)±a 



> 



e-? r np{S L n W) N x \k p \ [1 - (4 + if/fifi)]- 



The last inequality uses that |A«|(^ C + l) 2 //(/3) is the maximal number of sites in where it 
is not possible to insert a non-interacting critical droplet (recall (13. 12[) and note that a critical 
droplet fits inside an £ c x £ c square). According to Lemma [A. II in Appendix [S] we have 

MSinW)=M5x)[l + o(l)], (3-22) 

while conditions (|1.8fH.9j) imply that ^^{Sl) = fip(S)[l + o(l)]. Combining the latter with 
(|3.20ff3T2TT) . we obtain the desired lower bound. I 



3.2 Proof of Theorem fIT2l(b) 

We use the same technique as in Section \3.1\ which is why we only give a sketch of the proof. 

To estimate the average crossover time from 5^ C 5 to S C \C, we will use Proposition 12.11 
With A = S L and B = S C \C, (|2T0]) reads 

y 4 C \ C (a)E a (r s .\c) = CAHS l S c\ c) y ^(a)h SL ^ c (a). (3.23) 

The left-hand side is the quantity of interest in (11.171) . 

In Sections I3.2.1f437272l we estimate ^o-esue A i /3( <J )^5 L ,5 c \c( <T ) an d CAP(5l, S c \C). The 
estimates will show that 

r.h.s.([323D = ^-^e' 3r [l + (l)], (3 -» oo. (3.24) 
3.2.1 Estimate of £ ffe5uC W{°) h S L ,S-\c{°) 

Lemma 3.3 Y^aesyjc ^p( a ) h s L ,s-\c{^) = Vp(S)[l + o(l)] as (3 -> oo. 
Proof. Write, using (12. IB . 



H (a)h SL>S c\ e (a) = Hp{S L ) + y fi /3 (a)F a (T SL < t S c\ c ). (3.25) 
o-GSUC o-6(5\Sr)uC 

The last sum is bounded above by fj,p(S\Si,) + ^p{C). As before, fip{S\Si) = o{[ip{S)) as 
/? — ► oo. But (jl.35p and (|3.9|) imply that fJ,p(C) = o(fj,p(S)) as /? — > oo. I 
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3.2.2 Estimate of CAP(S L ,S C \C) 

Lemma 3.4 CAP(S,S C \C) = N 2 \Ap\e~ pr np(S)[l + o(l)] as (3 -> oo tfliifc iV 2 = f (24 - 1). 

Proof. The proof is similar as that of Lemma 13.21 except that it takes care of the transition 
probabilities away from the critical droplet. 

Upper bound: Recalling (|2.13fJ2~T5j) and noting that Glauber dynamics does not allow tran- 
sitions within C, we have, for all h: C — > [0, 1], 

CAP(<S L , S C \C) < CAP(S, S C \C) < £ HfiW) [ca(Ha) - l) 2 + c a (h(a) - 0) 2 ] , (3.26) 

o-ee 

where c CT = Y^neS c pk a -> V) an d ^ = J2 v es c \c c p( a > r ?)- The quadratic form in the right-hand 
side of (|3,26p achieves its minimum for h(a) = c a /(c a + c CT ), so 

CAP(S L , S C \C) < £C ff ^{a) (3.27) 



with C a = c a c CT /(c a + Co). We have 

(o») 



^8 



- • J_ V- e -/3^M 2 (I 4 + §(24 - 4)) ( 3 - 28 ) 
Z ? aev 



" N x 



where in the second line we use that C a = \ if a has a protuberance in a corner (2x4 choices) 
and C a = | otherwise (2 x (2£ c — 4) choices). 




Figure 10: Canonical order to break down a proto-critical droplet plus a double protuberance. In the 
first step, the double protuberance has probability g to be broken down in either of the two possible 
ways. The subsequent steps are deterministic as in Fig. [B] 

Lower bound : In analogy with (13.131) , denoting by P?* (x) the droplet obtained by adding a 
double protuberance at y along the longest side of the rectangle Rz c -i t i c (x), we define the set 
V L C S C \C by 

V L = {a\JPf y) {x): aeS L nW,x,yeAp,(x,y)±(r}. (3.29) 

As in (|3.15|) . we may choose any starting measure on T>l- We choose the flow as follows. For 
the first step we choose 

f(a', a) = \ u {a), a' £V L ,a£ C L , (3.30) 
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which reduces the double protuberance to a single protuberance (compare (|3.13p and (|3.29p ). 
For all subsequent steps we follow the deterministic paths j a used in Section 13.1.21 which 
start from 7 CT (0) = a. Note, however, that we get different values for the flows /(7 CT (0), 7cr(l)) 
depending on whether the protuberance sits in a corner or not. In the former case, it has only 
one possible antecedent, and so 



/(7.t(0),7.t(1)) = W*)> 
while in the latter case it has two antecedents, and so 

f(j a (0),j a (l)) = u (a). 



(3.31) 



(3.32) 



This time the terms k = and k = 1 are of the same order while, as in (|3.19p . all the 
subsequent steps give a contribution that is a factor 0(e _<5/3 ) smaller. Indeed, in analogy with 
(|3.17p we obtain, writing a ~ a' when cp{o\a) > 0, 

cap(s l ,s=\c) = cap(s'\c,s l ) > cap(v l ,s l ) 



/(a» + /(a,7,(l)) 



K-l 



WW) 



M°0 



f{ la {k), la {k + l)) 



cr'eV L <?ec L 



(3.33) 



[1 + (1)]^(C L ) 



[l + (l)]^(Ci) 



2t 



2L 



1 + 



1 4 1 

+ r 



\ 2 2 ^ c \ + 2 



JVi" 



Using (|3.2ip and the remarks following it, we get the desired lower bound. 



3.3 Proof of Theorem 

Write 

rtP M ^ ^s^V^ ^ 3 34 ^ 

= /i/3(5 L ) + E / x /3( fJ ) P ^( r 5L < TD M ). 

^ c m\Sl 

The last sum is bounded above by hp{S\Sl) + [ip{V c M \S). But ^p{S\Sl) = o([ip{S)) as 
/? — > oo by our choice of L in (|1.9p . while iip{V c M \S) = o(fj,p(S)) as /3 — ► oo because of the 
restriction i c < M2£ c — 1. Indeed, under that restriction the energy of a square droplet of size 
M is strictly larger than the energy of a critical droplet. 

Proof. The proof of Theorem II. 2l fc) follows along the same lines as that of Theorems II .2f a-b) 
in Sections EjJED The main point is to prove that CAP(S L ,V M ) = [1 + o(l)]CAP(S L , S C \C). 
Since CAP(Sl,T>m) < CAP(5l, S c \C), which was estimated in Section 13.21 we need only prove 
a lower bound on CAP(Sl,T>m)- This is done by using a flow that breaks down an M x M 
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droplet to a square or quasi-square droplet Ql in the canonical way, which takes M 2 — v(L) 
steps (recall Fig. and (|3.14|) ). The leading terms are still the proto-critical droplet with 
a single and a double protuberance. To each M x M droplet is associated a unique critical 
droplet, so that the pre- factor in the lower bound is the same as in the proof of Theorem ll.2l fb). 

Note that we can even allow M to grow with (3 as M = e°^>. Indeed, (|3.11ff3TT2|) show 
that there is room enough to add a droplet of size almost everywhere in A^, and the 
factor M 2 e~ & P replacing Ke~ 5p in (pO0|) still is o(l). I 



4 Proof of Theorem 11.41 

4.1 Proof of Theorem 11.4( a) 

4.1.1 Estimate of £ CT&S u(C\e+) M^s^s^uC+W 

Lemma 4.1 E cre su(C\C+) M cr ) /l s I ,,(«s c \e)uc+ ( a ) = Ppffli 1 + o(l)] as (3 -> oo. 
Proof. Write, using (|2.ip . 

/ i /3( (7 ) /l 5 i ,(5 c \C)UC+( cr ) 

CTe5 u(c\ C+ ) 

= H (S L )+ Yl ^W p »( T «t < r (5=\c)uc+)- 
CT e(5\5i,)u(c\c + ) 

The last sum is bounded above by hp(S\Sl) + ^ig(C\C + ). But ^p(S\Sl) = o(fj,p(S)) as 
/? — ► oo by our choice of L in (ll.32|) . In Lemma IB. 31 in Appendix IB. 31 we will show that 
Hp(C\C+) = o{fi p (S)) as p -> oo. 



4.1.2 Estimate of CAP(S L , (5 C \C) U C+) 

Lemma 4.2 CAP(5 L ,5 C \C) U C+) = N\Ap\ ^ e-^ fip(S)[l + o(l)] as [3 -» oo twtft iV 
§^-P 



Proof. The argument is in the same spirit as that in Section 13.1.21 However, a number of 
additional hurdles need to be taken that come from the conservative nature of Kawasaki 
dynamics. The proof proceeds via upper and lower bounds, and takes up quite a bit of space. 

Upper bound: The proof comes in 7 steps. 

1. Proto-critical droplet and free particle. Let C denote the set of configurations "in- 
terpolating" between C~ and C + , in the sense that the free particle is somewhere between 
the boundary of the proto-critical droplet and the boundary of the box of size Lp around the 
proto-critical droplet (see Fig. Ill|) . Then we have 

cap(5 l , (S C \C) U C+) < cap(5 U C- (S C \C) U C+) 

, min 2 W3{o)cp(a,o') [h(a) - h(a')] 2 . 

h: X P -.[0,1] 
h| SuC-= 1 ' ,l| (5 c \C)UC+= ^ 

(4.2) 
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Figure 11: Schematic picture of the sets S, C , C + defined in Definition 1 1 . 31 and the set C interpolating 
between C~~ and C + . 



Split the right-hand side into a contribution coming from a, a' G C and the rest: 



r.h.s.gj) = I + 7i(/3), 



where 



1= min \ Y] np{a)cp{a,a') [h(a) - h(a')f 

h: C— [0,11 



(4.3) 



(4.4) 



'ec 



and 7i (/?) is an error term that will be estimated in Step 7. This term will turn out to be small 
because iip{a)cp{a, a') is small when either a G X^ \c or a' G ^ n/3) \C. Next, partition C, 

C~ , C + into sets C(x), C~(x), C + (x), x G Ap, by requiring that the lower-left corner of the 
proto-critical droplet is in the center of the box B^^^^ix). Then, because cp(a, a') = when 
o G C(x) and a' G C{x') for some x ^ x' , we may write 



I A, 



mm 

h: C(0)-»[0,1] 



Vf3(<7)cf3{a,CT') [h(a) - h(a')} . 



'c-(o) = 



= l,h| c+(0) =0 a,a'GC(0) 



(4.5) 



2. Decomposition of configurations. Define (compare with (|3.6p ) 

C(0) = {al BLp ^ (0 y. cxGC(O)}, 

a G C(0)}. 



C(0) = {(71^ ^ (0)] c 



(4.6) 



Then every a G C(0) can be uniquely decomposed as a = a V<r for some a G C(0) and <r G C(0). 
Note that C(0) has = £ c (£ c — 1) + 2 particles and C(0) has np — K particles (and recall that, 
by the first half of (ll.35|) . rip — > oo as (3 — > oo). Define 



C fp (0) = {(J G C(0) : Hp(a) = Hp{a) + Hp{a)}, 



(4.7) 



i.e., the set of configurations consisting of a proto-critical droplet and a free particle inside 
-Bl^l^O) not interacting with the particles outside Sl^l^O). Write C fp '~(0) and C fp,+ (0) to 
denoting the subsets of C fp (0) where the free particle is at distance Lp, respectively, 2 from 
the proto-critical droplet. Split the right-hand side of (|4.5p into a contribution coming from 
a, a' G C fp (0) and the rest: 

r.h.s.(|43D = |A^| [II + 72^)], (4.8) 
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where 

11= min | V fxp(a)cp(a,a')[h(a)-h(a')] 2 (4.9) 

h: C f P(0)^[0,l] 

h l c fp,- ( o) sl >l C fP,+(o)-° -.-'^(0) 

and 72 (/3) is an error term that will be estimated in Step 6. This term will turn out to be 
small because of loss of entropy when the particle is at the boundary. 

3. Reduction to capacity of simple random walk. Estimate 

II = min 

h: C f P(0)->[0,l] 



h| cfp.-(o)= 1 '' l| ctp.+(o)= 



\ ^2 E M/3(£V a)cp(aV a, a' V a')[h(aV a) - h(a V a)} 2 



CTVCT,CT'Vo-'SC f P(0) 



< min 

g: <?(0)-*[0,l] 
9l C-(0) =1 ' 9, C+(0)=° 

\ £ A^ / 3('TV ( T)C /3 ( ( TV«T, ( 7'V ( 7)[ 5 («T)- 5 («T , )] 2 , 
o-eC(O) *,*'ec(0) : 

O"VlT,lT'VlTgC f P(0) 

where C (0), C(0) + denote the subsets of C(0) where the free particle is at distance Lp, 
respectively, 2 from the proto-critical droplet, and the inequality comes from substituting 

h(a Va) =g(&), a € 6(0), a € 6(0), (4.11) 

and afterwards replacing the double sum over cr, o 7 G C(0) by the single sum over a £ C(0) 
because 0/3(6" V a, a' V <r') > only if either <r = a' or it = a' (the dynamics updates one site 
at a time). Next, estimate 

r.h.s.glO]) 

^ E "L min I E e-^^^(a,a')b(a)-^)] 2 , 

^ 7W 9: C(0)-[0,1] T - ' 

C W C W CTV<T,<T'v«reC f P(0) 

(4.12) 

where we used Hp(a) = Hp(a) + Hp (a) from (|4.7p and write 0/3(6", 6"') to denote the transition 
rate associated with the Kawasaki dynamics restricted to 5x^^(0), which clearly equals 
cp(a V a, a' V 6") for every <t E C(0) such that a V 6", a' V <r G C fp (0) because there is no 
interaction between the particles inside and outside BLg,La(0)- The minimum in the r.h.s. of 
(|4.12p can be estimated from above by 

E v ^ cj ) ( 4 - 13 ) 

<xeP(o) 

with V(0) the set of proto-critical droplets with lower-left corner at 0, and 

Vp(a) = min \ £ [/(*) - /(x')] 2 , (4.14) 

/: Z2^[0,l] Z ^7^, 

/lp CT (o)-i./l[ Bi/3ii/3 (o)] C -o *'*'jf 



25 



where -Po-(O) is the support of the proto-critical droplet in a, and x ~ x' means that x and 
x' are neighboring sites. Indeed, (|4.13|) is obtained from the expression in (|4.12p by dropping 
the restriction a V a, a' V a G C fp (0), substituting 

g(P ff (0) U {x}) = f{x), a G V(0), x G B L ^(0)\P a (0), (4.15) 

and noting that c^(P cr (0) U {x},P a (0) U {x 1 }) = 1 when x ~ x' and zero otherwise. What 
(|4.13|) says is that 

Vp{a) = CAP(P CT (0), [B L0 . L0 (O)] C ) (4.16) 

is the capacity of simple random walk between the proto-critical droplet -Po-(O) in a and the 
exterior of Bl ,l„(O). Now, define 

Zf- K \<d)= Y, e- pHf,{&) . (4.17) 
o-eC(o) 

Then we obtain via (14.131) that 

r.h.B.CT < e^' g (wg ; £ V^a), (4.18) 
^8 CT eP(o) 

where T* = —U[(£ c — l) 2 + ^ c (^c — 1) + 1] is the binding energy of the proto-critical droplet 
(compare with (|1.33j) ) . 

4. Capacity estimate. For future reference we state the following estimate on capacities for 
simple random walk. 



Lemma 4.3 Let U C I? be any set such that {0} C U C Bi~^(0), with k G NU{0} independent 



o//3. Let 7 C Z 2 i»e any set swc/i that [B KL KL (0)] c C V C [B t L (0)] c , usi/i if G N 



independent of (3. Then 

CAP ({0}, [B KLptKLfl (0)] e ) < CAP (U, V) < CAP (5^(0), [B Lp , Lp (Q)] c ) . (4.19) 
Moreover, via (fQ9"HOU|) , 

CAP (S M (0), [^,^(0)] c ) = [1 + o(l)] lQg(gL/3) _ lQgfc = [1 + o(l)] ^, ^ - oo. 

(4.20) 

Proof. The inequalities in (|4.19p follow from standard monotonicity properties of capacities. 
The asymptotic estimate in (I4.20p for capacities of concentric boxes are standard (see e.g. 
Lawler |20j . Section 2.3), and also follow by comparison to Brownian motion. I 

We can apply Lemma 14.31 to estimate Vp(cr) in (|4.16p . since the proto-critical droplet with 
lower-left corner in fits inside the box -E>2£ c ,2£ c (0). This gives 

W°) = ^ [1 + o(l)], Va£ 7>(0), -+ oo. (4.21) 



Morover, from Bovier, den Hollander, and Nardi [7J, Lemmas 3.4.2-3.4.3, we know that iV 

M(3 



l'P(O)!, the number of shapes of the proto-critical droplet, equals N = \fi c {JP c — 1). 
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5. Equivalence of ensembles. According to Lemma iB.ll in Appendix [Bl we have 

1 ( / - = (pp) K M(S)[l + o(l)], P^oo. (4.22) 

This is an "equivalence of ensembles" property relating the probabilities to find np — K, 
respectively, np particles inside [B L/3i l (O)] c (recall (|4T6j) ). Combining (j4T2H43|) . (f4"3J) . P~8j) . 
(14~TU1) . ([4TT3D . (ITOD and (IT2TH021 .' we get 

CAP(S,C+) < 7l (/?) + |A /3 | 72 (/3)+iV|A /3 |— e-^ r ^(5)[l + o(l)], - oo, (4.23) 

where we use that T* + AK = T defined in (ll.33|) . This completes the proof of the upper 
bound, provided that the error terms 71 (/3) and 72(/3) are negligible. 

6. Second error term. To estimate the error term 72(/3), note that the configurations in 
C(0)\C fp (0) are those for which inside B L/3>L/3 (0) there is a proto-critical droplet whose lower- 
left corner is at 0, and a particle that is at the boundary and attached to some cluster outside 
BLf3,L (O)- Recalling (|4.5H4.9|) . we therefore have 

72(/3)< E M°)cp(<r,v')[h(a)-h(a')] 2 <6w(C(0)\C*(0)), (4.24) 

<rec(o)\e f p(o) o-'ec(o) 

where we use that h: C(0) — > [0, 1], pp(a)cp(a,a') = pp(cr) A pp(a'), and there are 6 possible 
transitions from C(0)\C fp (0) to C(0): 3 through a move by the particle at the boundary of 
BLf3,L (O) and 3 through a move by a particle in the cluster outside -Bl^l^O). Since 

Hp(a) > Hp(a) + H p (a) -U, a G C(0)\C fp (0), (4.25) 

it follows from the same argument as in Steps 3 and 5 that 

M/3 (C(0)\C fp (0)) < N e'^ (pp) K+1 n P {S) eF A(K - 1) [1 + o(l)], (4.26) 

where (pp) K+1 comes from the fact that rip — (K + l) particles are outside i^-i^-^O) (once 
more use Lemma IB. II in Appendix [B]), e@ u comes from the gap in (|4.25p . and A(K — 1) counts 
the maximal number of places at the boundary of B^^^^Q) where the particle can interact 
with particles outside Bl^^^Q) due to the constraint that defines S (recall Definition ll.3f) (a)). 
Since p e pu = o(l) by (fl~2?) . we therefore see that 72 (/3) indeed is small compared to the 
main term of (|4.23p . 

7. First error term. To estimate the error term 71 (/3), we define the sets of pairs of 
configurations 

1 1 = {(a, r,) G IX^} 2 -- aeS, V G S C \C}, 

t ^ ( 4 - 27 ) 
J 2 = {(<r, V ) G [X { p nf)) ] 2 : aeC, V e S C \C}, 

and estimate 

2 

7i(/?)<^E E ^(^)^(a ! r ? ) = iS(T 1 ) + iS(X 2 ). (4.28) 
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The sum can be written as 

E(2i) = |A/?|£ cp( V ,a)l{\suMv]nB L0 , L0 (O)\ = K} -^e-? H ^\ (4.29) 

where we use that pp(a)cp{a,rf) = fip(rj)cp(r],a), a, r] G and cp(rj,a) = 0, r? G 5 C \C, 

cr ^ P (recall Definition O^b)). We have 

fl/jfa) > fl^W) + Hp(r}) - kU, rj G S C \C, (4.30) 

where counts the number of pairs of particles interacting across the boundary of BLa,Lp{0)- 
Moreover, since rj ^ C, we have 

H p (fj)>T* + U. (4.31) 
Inserting (j43T)tiOTj) into (|439"j) . we obtain 

S(Jx) < \Ap\ n p (S) [1 + o(l)] Y^) K+k i< K ~ ^ k ~ 1)U (A ^ 

= |A /3 |e-^(5)[l + o(l)]e-^, 

where (pp) K+k comes from the fact that np — (K + A;) particles are outside Sx^_i,La-i(0) 
(once more use Lemma fB.ll in Appendix[B]), and the inequality again uses an argument similar 
as in Steps 3 and 5. Therefore £(Zi) is small compared to the main term of (|4.23p . The sum 
S(2"2) can be estimated as 

£(2" 2 ) = ^ Y Pp( a ) c p( a >V) 

a£Cr]eS c \C 

= \A p \ £ M/3 (<r) ^ c^,??) (4.33) 

treC(O) r;6<S c \C(0) 

< |A^| M/3 (C(0)) {e-^ + (4Lp) pp [1 + o(l)]}, 

where the first term comes from detaching a particle from the critical droplet and the second 
term from a extra particle entering Bl 0j £,^(O). The term between braces is o(l). Moreover, 
p,p(C(0)) = pp{C {p (0)) + pp{C{0)\C fp {0)). The second term was estimated in (jCTjl . the first 
term can again be estimated as in Steps 3 and 5: 

Z (n,3 ~ K) (0) 

Pp(C*(0)) = Y E HfiVV»)=Ne-Pr i _±l = Ne -pr^ {s)[1 + o{1)l 

aeC(0) *ec"(0) 
a-vo-6C f p(o) 

(4.34) 

Therefore also S(X 2 ) is small compared to the main term of (|4.23l) . 

Lower bound : The proof of the lower bound follows the same line of argument as for Glauber 
dynamics in that it relies on the construction of a suitable unit flow. This flow will, however, 
be considerably more difficult. In particular, we will no longer be able to get away with 
choosing a deterministic flow, and the full power of the Berman-Konsowa variational principle 
has to be brought to bear. The proof comes in 5 steps. 

For future reference we state the following property of the harmonic function for simple 
random walk on Z 2 . 
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Lemma 4.4 Let g be the harmonic function of simple random walk on B2L g ,2L g (0) (which is 
equal to 1 on {0} and on [B2L B ,2L g (0)] c )- Then there exists a constant C < oo such that 

J2[g(z)-g(z + e))+ < C/L p Vz G [B L ^ Lp (0)] c . (4.35) 

e 

Proof. See e.g. Lawler, Schramm and Werner [21j . Lemma 5.1. The proof can be given via the 
estimates in Lawler [20], Section 1.7, or via a coupling argument. I 



1. Starting configurations. We start our flow on a subset of the configurations in C + that 
is sufficiently large and sufficiently convenient. Let C C + denote the set of configurations 
having a proto-critical with lower- left corner at some site x G A^, a free particle at distance 2 
from this proto-critical droplet, no other particles in the box B2t, B ,2L B { x )i an d satisfying the 
constraints in Sl, i.e., all other boxes of size 2L ( g carry no more particles than there are in a 
proto-critical droplet. This is the same as C + , except that the box around the proto-critical 
droplet has size 2Lp rather than La. 

Let K = i c (i c — 1) + 2 be the volume of the critical droplet, and let S^ 13 K ^ be the 
analogue of S when the total number of particles is ng — K and the boxes in which we count 
particles have size 2Lp (compare with Definition II. 3D . Similarly as in A3. 1TH . our task is to 
derive a lower bound for CAP(S L , (5 C \C)UC+) = CAP((S C \C) UC + ,S L ) > CAP(C L ,S L ), where 
C L CC^ CC + defined by 

C L = {aU P [y) (x, z): a£ S { ^' K \ x, y G Ag, (x, y, z)±a} (4.36) 

is the analog of (|3,13p . namely, the set of configurations obtained from S% K ^ by adding a 
critical droplet somewhere in A^ (lower-left corner at x, protuberance at y, free particle at z) 
such that it does not interact with the particles in a and has an empty box of size 2Lg around 
it. Note that the np — K particles can block at most np(2Lp) 2 = odA^I) sites from being the 
center of an empty box of size 2L ( g, and so the critical particle can be added at [A^l — odA^I) 
locations. 

We partition Cl into sets Cl(x), x G A^, according to the location of the proto-critical 
droplet. It suffices to consider the case where the critical droplet is added at x = 0, because 
the union over x trivially produces a factor |A^|. 

2. Overall strategy. Starting from a configuration in Cl(0), we will successively pick K — L 
particles from the critical droplet (starting with the free particle at z at distance 2) and 
move them out of the box 5^,^(0), placing them essentially uniformly in the annulus 
-B2L /3 ,2L (3 (0)\-Bl^,l /3 (0). Once this has been achieved, the configuration is in Sl. Each such 
move will produce an entropy of order L^, which will be enough to compensate for the loss 
of energy in tearing down the droplet (recall Fig. U]). The order in which the particles are 
removed follows the canonical order employed in the lower bound for Glauber dynamics (recall 
Fig. [6]). As for Glauber, we will use Proposition 12.41 to estimate 



CkV{C Ll S L )>\A p \ E p/ W 



t( 7 ) 

E 

o-eC L (0)7: 70=c k=0 



/(7fc,7fe+i) 



1 -1 



^/3(7fc) c /3(7fc ) 7fc- 



(4.37) 



for a suitably constructed flow / and associated path measure W, starting from some initial 
distribution on Cl(0) (which as for Glauber will be irrelevant), and t(j) the time at which 
the last of the K — L particles exits the box Bl„^A<S) . 
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The difference between Glauber and Kawasaki is that, while in Glauber the droplet can be 
torn down via single spin- flips, in Kawasaki after we have detached a particle from the droplet 
we need to move it out of the box -6^,1^(0), which takes a large number of steps. Thus, ri^f) 
is the sum of K — L stopping times, each except the first of which is a sum of two stopping 
times itself, one to detach the particle and one to move it out of the box 5/^^(0). With each 
motion of a single particle we need to gain an entropy factor of order close to l/pp. This will 
be done by constructing a flow that involves only the motion of this single particle, based on 
the harmonic function of the simple random walk in the box -E^L^, 2^(0) up to the boundary 
of the box Bl 0) l„(O). Outside BLg,Lp(fy the flow becomes more complex: we modify it in 
such a way that a small fraction of the flow, of order L^ 1+e for some e > small enough, is 
going into the direction of removing the next particle from the droplet. The reason for this 
choice is that we want to make sure that the flow becomes sufficiently small, of order L^ 2+e , 
so that this can compensate for the fact that the Gibbs weight in the denominator of the lower 
bound in (I2.20p is reduced by a factor e~^ u when the protuberance is detached. The reason 
for the extra e is that we want to make sure that, along most of the paths, the protuberance 
is detached before the first particle leaves the box -B2L /3 ,2L (3 (0). 

Once the protuberance detaches itself from the proto-critical, the first particle stops and 
the second particle moves in the same way as the first particle did when it moved away 
from the proto-critical droplet, and so on. This is repeated until no more than L particles 
remain in 5^,^(0), by which time we have reached Sl- As we will see, the only significant 
contribution to the lower bound comes from the motion of the first particle (as for Glauber), 
and this coincides with the upper bound established earlier. The details of the construction 
are to some extent arbitrary and there are many other choices imaginable. 

3. First particle. We first construct the flow that moves the particle at distance 2 from 
the proto-critical droplet to the boundary of the box jBl^i^O). This flow will consist of 
independent flows for each fixed shape and location of the critical droplet. This first part of 
the flow will be seen to produce the essential contribution to the lower bound. 

We label the configurations in Cl(0) by a, describing the shape of the critical droplet, as 
well as the configuration outside the box BzLg^Lpify, and we label the position of the free 
particle in a by Z\{a). 

Let g be the harmonic function for simple random walk with boundary conditions on 
[-S2L (3 ,2L /3 (0)] C and 1 on the critical droplet. Then we choose our flow to be 

f(a(z), a(z')) = { Cl [9{Z) ~ 9{Z + if Z ' = Z + 6 ' l|e|1 = h (4.38) 

I 0, otherwise, 

where cr(z) is the configuration obtained from a by placing the first particle at site z. The 
constant C\ is chosen to ensure that / defines a unit flow in the sense of Definition 12.31 i.e., 

Y, C± J2{9(zi(a))-g(z 1 (a) + e)]=C 1 £ CAP (P ff (0), [B 2L ^ 2L ^)f) = 1, (4.39) 

aeC L (0) zi(<r),e <reC L (0) 

where -Po-(O) denotes the support of the proto-critical droplet in a, and the capacity refers to 
the simple random walk. 

Now, let z 1 (k) be the location of the first particle at time k, and 

r 1 = inf{A: G N: z\k) € [B L ^(0)] C } (4.40) 
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be the first time when, under the Markov chain associated to the flow /, it exits -Bl^l^O). 
Let 7 be a path of this Markov chain. Then, by (|4.38fj4T3"9"jh we have 



E 

fc=0 



/(7fc,7fc+i) 



^/3(7fc)c/3(7fc,7fcH 



^(7o) 



(4.41) 



where the sum over the g's is telescoping because only paths along which the (^-function 
decreases carry positive probability, and 0,3(7^,7^+1) = 1 for all < k < r 1 because the first 
particle is free. We have (/(^(O)) = 1, while, by Lemma 14.41 there exists a C < 00 such that 



Therefore 



E 

k=0 



g{x) <C/ log Lp, x € [B LfhLf ,(0)] 
1 

/(7fe)7fc+i) Ci 



M/3(7fc)c/3(7fc,7fc+i) ^(70) 



l + o(l)]. 



Next, by Lemma 14.31 we have 



4tt 



CAP(P CT (0), [B 2Lf3>2L0 (O)) c ) = ^[l+o(l)], rredH)). > x. 



(4.42) 
(4.43) 

(4.44) 



(because {0} C P ff (0) C B 2 i> c . 2l > c (0) for all cr G C L (0)). Since TV = |C L (0)|, it follows from 
<|Q5|> that 

k = N jk [l+o{l)l (445) 

and so (|4.43j) becomes 



E 

fc=0 



/(7Jfc>7fc+i) 



M/3(7fc)c/3(7fc,7fc+i) 
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W(-Yo)N— [l + o(l)]. 



/3 — > OO, 



(4.46) 



This is the contribution we want, because when we sum (|4.46|) over 70 = a G Cl(0) (recall 
(I4.37P ). we get a factor 

^(C L (0)) = e"^ ^(S) [1 + o(l)]. (4.47) 

To see why (I4.47D is true, recall from (14.360 that Cl(0) is obtained from S^ 13 K ^ by adding 
a critical droplet with lower- left corner at the origin that does not interact with the rip — K 
particles elsewhere in A^. Hence 



MCl(O)) 



-/3r* z£ 



7 ( n p) 



(4.48) 



where zf K) (0) is the analog of K '(0) (defined in <^T7j> ) obtained by requiring that 
the rip — K particles are in [Re c ,e c (0)] c instead of [Bi„^„(0)] c . However, it will follow from 
the proofs of Lemmas IB . H4b721 in Appendix [B] that, just as in (|4.22|) . 



{ P p) K Vp(S) [l + o(l)], /3roo, 



(4.49) 
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which yields (|4.47p because r = T* + K A. For the remaining part of the construction of the 
flow it therefore suffices to ensure that the sum beyond r 1 gives a smaller contribution. 

4. Second particle. Once the first particle (i.e., the free particle) has left the box Sl^l^O), 
we need to allow the second particle (i.e., the protuberance) to detach itself from the proto- 
critical droplet and to move out of Bl /3 ^ /3 (0) as well. The problem is that detaching the 
second particle reduces the Gibbs weight appearing in the denominator by e~ u @, while the 
increments of the flow are reduced only to about 1/Lp. Thus, we cannot immediately detach 
the second particle. Instead, we do this with probability L^ 1+e only. 

The idea is that, once the first particle is outside £?£„£„(()), we leak some of the flow that 
drives the motion of the first particle into a flow that detaches the second particle. To do 
this, we have to first construct a leaky flow in -B2L (3 ,2L^(0)\-Bl /3i l (3 (0) for simple random walk. 
This goes as follows. 

Let p(z, z + e) denote the transition probabilities of simple random walk driven by the 
harmonic function g on -E^L^L^ (0). Put 

P[ ' + j "\ (1- Lp 1+e )p(z,z + e), ifze B 2L0 , 2Lfj (O)\B L „L p (.O). ^ 

Use the transition probabilities p(z, z + e) to define a path measure P. This path measure 
describes simple random walk driven by g, but with a killing probability L^ 1+e inside the 
annulus B 2 l ,2L^)\Bl ,l p {^)- Put 

k(z, z + e)=Y, £(7)l(*,*+e)ey, z G #2^,2^(0). (4.51) 

7 

This edge function satisfies the following equations: 

• k(z, z + e) = [g(z) - g(z + e)]+, 

if z G B L ^ LfJ (0), 

• k(z, z + e) = 0, 

if z G j B 2 l,,2L,(0)\ j Bl /3 ,l /3 (0) and [g{z) - g(z + e)]+ = 0, (4.52) 

• (f - L p 1+e ) ^2 H z + e > z)l g ^ z+e) ^ g{z)> Q = K z , z + e)t giz) _ g{z+e)>0 

e e 

if z e B 2L ^ 2Lg (0)\B L ^ Lg (0). 

Note that inside the annulus B 2 Ln,2La (0)\5x^,Lg(0) at each site the flow out is less than the 
flow in by a leaking factor 1 — L7, 1+e . We pick e > so small that 



e^ u is exponentially smaller in f3 than L 2 g e , (4.53) 

(which is possible by (jl.27p and (jl.291fl~30|) ). The important fact for us is that this leaky flow 
is dominated by the harmonic flow associated with g, in particular, the flow in satisfies 

Y J k(z + e,z)<Y[g(z + e)-g(z)} + VzG%,, 2 ^(0), (4.54) 

e e 

(and the same applies for the flow out). This inequality holds because g satisfies the same 
equations as in (|4.5UH4"3T]) but without the leaking factor 1 — L^ 1+e . 
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Using this leaky flow, we can now construct a flow involving the first two particles, as 
follows: 

• f(a(zi,a),a(z 1 + e,a)) = Cik(z 1 ,zi + e), 

if z\ G B 2 l p ,2L (O), 

• f(a(z 1 ,a),a(z 1 ,b)) = CiL^ 1+e ^ k(zi, z\ + e), 

e 

if *1 e52L s ,2L fl (0)\flL s ,L fl (0), (455) 



• /(o-(^i,^2),cr(2i,z 2 +e)) = < CiL. i+e V + e) }> \g{z 2 ) - g(z 2 + e)]+, 



^2k(zi,zi + e) j 



if *i E B 2L ^ 2Lf3 (O)\B L0)L0 (O),z 2 G B L0}L0 (O)\P a (O). 

Here, we write a and b for the locations of the second particle prior and after it detaches itself 
from the proto-critical droplet, and a(z\, z 2 ) for the configuration obtained from a by placing 
the first particle (that was at distance 2 from the proto-critical droplet) at site z\ and the 
second particle (that was the protuberance) at site z 2 . The flow for other motions is zero, and 
the constant C\ is the same as in (|4.38fPO§j) 

We next define two further stopping times, namely, 

C 2 = inf{A: G N: z 2 ( 7fc ) = b}, (4.56) 

i.e., the first time the second particle (the protuberance) detaches itself from the proto-critical 
droplet, and 

r 2 = inf{A; G N: z 2 ( lk ) G [Bl,,l,(0)] c }, (4.57) 

i.e., the first time the second particle exits the box Bl^^^^O). Note that, since we choose the 
leaking probability to be L~ 1+e , the probability that Q 2 is larger than the first time the first 
particle exits -B2L /3 ,2L /3 (0) is of order exp[— L e g] and hence is negligible. We will disregard the 
contributions of such paths in the lower bound. Paths with this property will be called good. 

We will next show that (|4.4ip also holds if we extend the sum along any path of positive 
probability up to £ 2 . The reason for this lies in Lemma flow- lb. 11. Let 7 be a path that has 
a positive probability under the path measure W associated with / stopped at r 2 . We will 
assume that this path is good in the sense described above. To that end we decompose 



f(jk,lk+l) 



_ = y- /(7fc,7fc+i) /(7fc»7fc+i) /Q7fc 1 7fc+i) < ^ 4 " 58 ' ) 
^ l^p{lk)cp{-ik,lk+i) k £^ +1 ^/3(7fc)c/3(7fc,7fc+i) fe 4^-i Wil^cpi-ykilk+i) 

= 1 + 11 + III. 

The term / was already estimated in (14,411 [4~4"T1) . To estimate II, we use (I4,42p and (I4,54fj4~55j) 
to bound (compare with (|4.4ip ) 

, ;<Ci jMj)< C ,M, (4.59) 
M/H7o) Hp (70) 
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which is negligible compared to I due to the factor C/logL^. It remains to estimate III. 
Note that 



_2 



/i/3(7c 2 -i)c/3(7c 2 -i>7c 2 ) fc 4^ 2 M/3(7fc)c/3(7fe>7fc+i)' 

The first term corresponds to the move when the protuberance detaches itself from the proto- 
critical droplet. Its numerator is given by f(a(zi,a),a(zi,b)) (for some z\ 6 [-Bz^.Z^ (0)] c ) 
which, by Lemma [Hand (f4T54H435l) . is smaller than dL^'CL' 1 = dCL p 2+e . On the 
other hand, its denominator is given by 

zi(7 C 2_ 1 )c /3 (7 C 2_ 1 ,7 C2 ) = zi /3 (7o)e- c//3 . (4.61) 

The same holds for the denominators in all the other terms in III, while the numerators in 
these terms satisfy the bound 

/(7*,7*+i) < CrC L- 2 ^[g(z 2 { lk )) - g{z 2 { lk+l ))]. (4.62) 

Adding up the various terms, we get that 

III < -^V^(l + [g(z 2 (C 2 )) - g (z 2 {r 2 )}) < J^Lf^fF. (4.63) 
At/3 (7o ) p V J np(no) p 



The right-hand side is smaller than I by a factor L /3 2+e e^^, which, by (|4.53p . is exponentially 
small in (3. 

5. Remaining particles. The lesson from the previous steps is that we can construct a 
flow with the property that each time we remove a particle from the droplet we gain a factor 
Lp 2+e , i.e., almost e _A ^. (This entropy gain corresponds to the gain from the magnetic field 
in Glauber dynamics, or from the activity in Kawasaki dynamics on a finite open box.) We can 
continue our flow by tearing down the critical droplet in the same order as we did for Glauber 
dynamics. Each removal corresponds to a flow that is built in the same way as described in 
Step 4 for the second particle. There will be some minor modifications involving a negligible 
fraction of paths where a particle hits a particle that was moved out earlier, but this is of no 
consequence. As a result of the construction, the sums along the remainders of these paths 
will give only negligible contributions. 

Thus, we have shown that the lower bound coincides, up to a factor 1 + o(l), with the 
upper bound and the lemma is proven. I 



4.2 Proof of Theorem EKb) 

The same observation holds as in (|3.34p . 

Proof. The proof of Theorem ll.4f b) follows along the same lines as that of Theorem 1 1.4 I f a). 
The main point is to prove that CAP(£> A /,Sl) = [l + o(l)]CAP(C + ,S L ). Since CAP(S L ,V M ) < 
CAP(5z,C + ), we need only prove a lower bound on CAP(T>m, Sl)- This is done in almost 
exactly the same way as for Glauber, by using the construction given there and substituting 
each Glauber move by a flow involving the motion of just two particles. 

Note that, as long as M = e°^\ an M x M droplet can be added at |A^| — o(|Ag|) locations 
to a configuration a G S (compare with (14.360 ). The only novelty is that we have to eventually 
remove the cloud of particles that is produced in the annulus B2l ,2Lb(O)\Bl 0j l„(O). This is 
done in much the same way as before. As long as only particles have to be removed, 
potential collisions between particles can be ignored as they are sufficiently unlikely. I 
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A Appendix: sparseness of subcritical droplets 

Recall Definition II. If a) and (|3.11I[3TT2j) , In this section we prove (|3,22p . 
Lemma A.l lim^ J log ^1 = -oo. 

Proof. We will prove that lim^oo jflog (Jtp(S\W) / (ip(B) = — oo. Since B6 5, this will prove 
the claim. 

Let f{(5) be the function satisfying (|3.11|) . We begin by noting that 

Hp(S\W) < n fi {T) with T= {a G S: |supp[C B (a)]| > |A^|//(/3)}, (A.l) 

because the bootstrap percolation map increases the number of (+l)-spins. Let T>{k) denote 
the set of configurations whose support consists on k non-interacting subcritical rectangles. 
Put C\ = (£ c + 2)(£ c + 1). Since the union of a subcritical rectangle and its exterior boundary 
has at most C\ sites, it follows that in X there are at least \Ap\/Cif((3) non-interacting 
rectangles. Thus, we have 

MX)< £ F(k) with F(k) = — e-? H ^\ (A.2) 

CiUfi} c( CT )e-D(fc) 

where K max <\A p \. 
Next, note that 

F(k) < (2 Cl ) k -J- V e -^W. (A.3) 

Since the bootstrap percolation map is downhill, the energy of a subcritical rectangle is 
bounded below by C2 = 2 J — h (recall Fig. [9]), and the number of ways to place k rect- 
angles in A^ is at most ('^f'), it follows that 



Pk], 



F(k) < 2 Clfc ^(B)e- C ^ k < 2 c ^(C ie f(P)) k ^(B)e- c ^ k < fi (B) exp[-±C : 

(A.4) 

where the second inequality uses that k\ > k k e k , k G N, and the third inequality uses that 
f(p) = e °^\ We thus have 



14 1 i^JA^I 



F(k) <2 / u /3 (B)/(/3)M exp 



(A.5) 



which is the desired estimate because f((3) tends to infinity. I 

B Appendix: equivalence of ensembles and typicality of holes 

For m G N, let 

SM = {a G Xjp : |suppH n B L ^{x)\ < 4(4 - 1) + 1 Vx G A^} (B.l) 
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and 

c (m) (fl) = K v .(») :ffe5W }' 

zj™\o)= Yl e ~ pH{a) - (B ' 2) 

ff £CH(o) 

The latter is the partition sum restricted to -Bl^l^O) when it carries m particles. In Ap- 
pendix [Bj] we prove a lemma about ratios of partition sums that was used in (|4.22p . (|4.26p . 
<|Q2|> and fliagl) . In Appendix E2 we prove that lim i g^ oo /x i g(5(0))//i i g(5) = 1, which is 
needed in the proof of this lemma. 

B.l Equivalence of ensembles 
Recall (fQ2l . (IPD and ([OTP . 

Lemma B.l Z^ 8 ~ s \o)/Z^ 8) = {pp) s pp{S) [1 + o(l)] as /3 -> oo /or a// s G N. 

Proof. The proof proceeds via upper and lower bounds. 
Upper bound: Let 

5(0) = {a G 5 : suppH fl B^,^ (0) = 0} . (B.3) 

Write 

«*4e e (";) ^"' (svc) . (b.4) 

Z /3 d-eC(O) Cc[-B i/3ii/3 (0)]c\su P p[*] v ' 
Kl=« 

This relation simply says that rip particles can be placed outside BLp,L«(P) by first placing 
rip — s particles and then placing another s particles. Because the interaction is attractive, 
we have 

Hp(&vQ < Hp(a) + Hp(() and H p (Q < 0, V<t,C- (B.5) 

Consequently, 

^(o))>(7)"-^r E o-""'" E w<*»- 

Z /3 <t6C(0) Cc[s i(3iL(3 (0)] = \su PP [*] 

We next estimate the second sum, uniformly in a. When we add the s particles, we must 
make sure not to violate the requirement that all boxes Bl^^Jx), x G A^, carry not more 
than K particles (note that 5(0) C S and recall Definition 11.3( a)). Partition Ap\Br, 0t i g (O) 
into boxes of size Lp. The total number of boxes containing K particles is at most rip/K. 
Hence, the total number of sites where we cannot place a particle is at most (np/K)(3Lp) 2 . 
Therefore 



CC[B L/3 , L/3 (0)]c\{a} 

ICI=« 



V- -„ (\Kp\B L ^)\-np-(np/K){3Lpf\ 

2^ %vc&S(0)} > [ s J ' V it. (B.7) 



36 



But ji/jLj = o{np/p p ) = o{\kp\) and L 2 = o(l/p p ) = o(\Ap\) by (fL22l) and OSIOO]) . and 
so the right-hand side of (|B.7p equals [1 + o(l)] |A^| s /s! as (3 — > oo. Since the binomial in 
(IBTBT) equals [1 + o(l)] (np) a /a\ with = \pp\Ap\], we end up with (recall (|4.17p ) 

^(5(0)) > g (nfl) V ( P/} )- [1 + (1)], (B.8) 

or 

g ( } < (W) s ^(5(0)) [1 + o(l)]. (B.9) 

Z /3 

Since 5(0) C 5, this gives the desired upper bound. 

Lower bound: Return to (jB.4|) . Among the s particles that are added to [i?£, /5j L /3 (0)] c , let si 
count the number that interact with the rag — s particles already present and S2 the number 
that interact among themselves, where s\ + S2 < s. We can then estimate 



Z /5 *ec(0) V 7 o S .i+? a < 



e l{|Cn<9o-|=si} I{s2 interacting particles in ^{ctVCGS(O)} 



CC[S L/3ii/ j(0)] c \ S upp[«T] 



Z f3 



e «-**<•> e 

x ^ ] e ^ ''^ l{|Cn9o-|=si} l{s 2 interacting particles in C} l{o-VC6<S(0)}' 



CC[S Lj3ji)3 (0)] c \supp[a] 
ICI=« 

(B.10) 

where in the second inequality we estimate the term with si = S2 = by using the result for 
the upper bound. We will show that the other terms are exponentially small. 

For fixed a, we may estimate the last sum in (jB.lOp as 

^ ^ e ^ ^ l{|(n96-|=si} l{s 2 interacting particles in ^{<5-VCe<S(0)} 



CCI-B^ , (0)] c \supp[S] 



(B.ll) 

< I A I s — -si — s 2 (An \ s ~l \ ^ p —PH( a ) 1 

— l Ji pl \^ n /3) / j t- -"-{^2 interacting particles in <r} - 

Indeed, |A / g| s_,Sl_ ' S2 bounds the number of ways to place s — s% — S2 non-interacting particles, 
and (4n ( g) Sl the number of ways to place s% particles that are interacting with the rip — s 
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particles already present. Next, we write 

E p -/3H(a) t 
{«2 interacting particles in a] 

cr e 5< s 2) 

= EE E E e "' 3E ' =lfl(C ' ) . 

m=l jf=l 2<fei,...,fc i <K c=u 3 =iC . 

which is a cluster expansion of the partition function (with non-interacting clusters C{, all 
of which have size < K = £ c (i c + 1) + 1)- By a standard isoperimetric inequality we have 
H{Ci) > H^., with the latter denoting the energy of a droplet of ki = \Ci\ particles that is 
closest to a square or quasi-square. Hence 

I a |-s 2 V -PH(a) -« 

M^/Jl / e {s2 interacting particles in cr} 

S2 m / \ 

<^r>EE E E 1 

m=l j=l 2<fci fcj<Jf \ 0=u j 0j / 

Ef =1 fej=s2 icy^vi 
S2 m 

< C\A?\- S2 EE E ' °" IAb| j ' 

2<fe 1 ,...,*y <K 

fe i= s 2 



(B.13) 



-fi E^iEHii+C**-!)^- 1 log ia^h 



,*y<K 



S2 m 



-EE E e- /3E -i^ +(fcl " 1)A] , 

m=l 3=1 2<k 1 ,...,k j <K 

Ef =1 *>i=*2 

where in the last inequality we insert the bound log |A#| > A, which is a immediate from 
(TL22D and (05]) . 

Now, +k{A is the grand- canonical energy of a square or quasi-square with ki particles. 
It was shown in the proof of Proposition 2.4.2 in Bovier, den Hollander and Nardi [7J that 
this energy is > U \fk{ for 1 < k% < i.e., for a droplet twice the size of the proto-critical 
droplet. Since 2U > A, we therefore have that H/-. + (ki — 1)A > when ki > 4. Since A > U, 
H2 = —U and H3 = —2U, we have that also the terms with ki = 2,3 are > 0. Consequently, 
there exist e > and a constant C that is independent of (3 such that 

I A |- s 2 p ~PH(a) 1 -/3H(cr) < r -pt (VI 1 A} 

r 4 -^! / y e {«2 interacting particles in cr}'' _^ L/ c . ^-D.l'lj 

cr e 5< a 2) 

Combining (|B.10l4T37i"TT) and (jB. 14[) . we see that the correction term in (|B,10p is 



^(«S(0))-[1 + O (l)] P (n/)) (ft0~ 

y{ n f3- s ), 



Z v "p J (o) 

<c[i+ (i)] \ p p y s e (^Vr 1 ^. 



(B.15) 
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Since A > U, we have e u ^pp < 1 and so the sum is o(l). Hence 

& {np) > (Pfi)' ^(5(0)) [1 + (1)]. (B.16) 

Z /3 

The claim now follows by using Lemma lB.21 below. I 

B.2 Typicality of holes 

Lemma B.2 lim^oo ^g(5(0))//x^(5) = 1. 

Proof. Since 5(0) C S, we have ^(5(0)) < fJ-p(S). It therefore remains to prove the lower 
bound. Write 

M/3 («S) = ^(S(0)) 

K -/3H(t)V0 
+ E E E - ~ l{ S upp[T,]cB L0>Lf} (0)}l{supp[C]C[B £ ^, ifl (0)] c } 

<^(5(0))+7i(/3)+72(/3), 

(B.17) 

where 

e -/3[H( v )+H(Q] 

7l(/?)=2^ 2^ 2^ 1 {supp[r,]cB L;3 , L/3 (0)}l{supp[C]C[B L/3 , L/3 (Or} (B.18) 

Tn=1 ue*j m) ce^ n - m) Z /3 

and 72 (/?) is a term that arises from particles interacting accross the boundary of Bl 0! i,^(O). 
We will show that both 71 and 72 (/3) are negligible. 

Estimate, with the help of (|B.9j) (and recalling (TB.lHB.2j) ). 

1 V" P -PH{ri) 



m=l ^3 ^e^f" 1 ' 

K 

[1 + O (l)]^(5(0))^(^) m 2 1 {su PP [ J3 ]cB £ „ i ,(o)} (B>19) 

m=l r;6<S( m ) 
A' m 

[i+oWi^iE^rE E E e-^^, 



m=l 3=1 2<ft!,...,fcj<K c= j c 
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where the last equality is a cluster expansion as in (|B.12p . Using once more the isoperimetric 
inequality, we get (recall ([1,291) ) 



K 



7i(/?)<[l + O (l)]M<5(0)) E^ m E E e~^ =1 ^) £ 1 

m=l j=l 2<k 1 ,...,k j <K \ c=u i =iC . 



< cr^(0))X; for E 



i=i 2<fci,...,fc f <K (B.20) 



if m 

= C^(«S(0)) £ E E e-^LJ^+^A-CA-^)] 

m=l j = l 2<*i,...,fcj<K 

< Ciif,(S(0)) e"* 

for some e > and constants C,C < oo that are independent of /?. 
Estimate, with the help of (|B.9|) . 

K m y{rif3-(m+k)) 

7 2 (/3)<E E ^ W E^Wy f iV- 
m=l ijeS( m ) fc=1 /3 



<E E e-WJe^l^^^^r' M^)) [l + o(l)] 
m=l n&s^ m ) fc=i 

if m 

< [i+ (i)i^(5(o))^w m E e "^ (J7) E e " /3fc(A ^ )l wpwc^ /3 ,^(o)}> 

m=l r)£S( m ) fc=l 

(B.21) 

and we can proceed as (IB,19flB~20]l to show that this term is negligible. I 

B.3 Atypicality of critical droplets 

The following lemma was used in Section 14,1.11 
Lemma B.3 lim j 8_ +00 ^p(C\C + ) / ^p(S) = 0. 
Proof. Similarly as in (1B.17|1 . we first write 

e ~P[H(r,)+H(C)} 

= |A^| 7(/3) + IA/jI 2^ 2^ l{supp[r,]cB i(3 , L(3 (0)}l{supp[C]C[B i/3 , L(3 (Or}- 

(B.22) 

with 7(/3) a negligible error term that arises from particles interacting accross the boundary 
of B Lj3}L0 (O). We then proceed as in ijElsHB^Ojl . obtaining (r = T* + KA) 



r.h.s.dESI < N \A p \ {p p ) K ^(5(0)) [1 + o{l)\ 



N\A P \ e-^ ^(S)[l + o(l)}, (3 



(B.23) 



oo, 
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which is o(fj,p(S)) by ()1.35|) . 



C Appendix: Typicality of starting configurations 

In Sections IC.lHC.2l we prove the claims made in the remarks below (jl.9jl . respectively, (|1.32j) . 



C.l Glauber 

Proof. Split 



s = s L u(s\s L ) = s L uu >L , (c.i) 



where U > l C S are those configurations a for which Cb(<t) has at least one rectangle that is 
larger than Ql(0). We have 

C B {a)= |J # 4(a y 2(a;) 0), (C.2) 

x6X(<r) 

where X(cr) is the set of lower-left corners of the rectangles in Cb(ct), which in turn can be 
split as 

X{a) = X >L (a)UX^ L (a), (C.3) 

where X >L (a) labels the rectangles that are larger than Ql(0) and X—(a) labels the rest. 
Let g\a denote the restriction of a to the set A C Z 2 . Then, for any x E X(a), we have 

F W = F ( ff k(-)W-)) +F H^wW-))' (c - 4) 

because the rectangles in Cb{o~) are non-interacting. Since for cr G U > l there is at least one 
rectangle with lower-left corner in X >L (a), we have 

H(3{U>l)< E E 1 {xex>L(a)} 



< e -^+i ^ E V v '••■•17:' ' 1 



1 -phL\ r c (x) ) 

1 {ie^ >i W} ; 



(C.5) 

where r^+i is the energy of Ql + i(0). In the last step we use the fact that the bootstrap map 
is downhill and that the energy of Ql(0) is increasing with L. Since the energy of a subcritical 
rectangle is non-negative, we get 

M/9(W>l) < N L+1 e-^+i \Ap\ n (S) (C.6) 

with Nl + \ counting the number of configurations with support in Q^ + i(0). 

On the other hand, by considering only those configurations in U > l that have a Ql+i(0) 
droplet, we get 

^(U>l) > N L+1 e-^ |A^| M ^ +l(0)]C (5), (C.7) 
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where the last factor is the Gibbs weight of the configurations in S with support outside 
[Ql+i(0)] c . It easy to show that nf L+im \s) = np(S)[l + o(l)] as (3 -» oo and so 

fip{U >L ) > N L+1 e~^^ \Ap\fxp(S)[l + o(l)], /3^oo. (C.8) 
Combining (jC.6HC.7|) . we conclude that lim^oo [ip(U > L) / np(S) = if and only if 

lim |A«| e- TL+1 = 0. (C.9) 



C.2 Kawasaki 

Proof. Split 

s = s L u(s\s L ) = s L uu >L , (c.io) 

where U>l C S are those configurations a for which there exists an x such that |supp [a] n 
S L/3 , L/3 (x)| >L. Then 

< E E E W( a ) 1 {\™pp{*}nB L L (x)\= m } = \^\[<p{f3)+j(P)], (C.ll) 



x£Af) <tGS m=L+l 



where 



e -|9[ff(,)+H(()] 

<P(P)= E E E Z(^) l{suppMcS i/3 , i/3 (0)} l{supp[C]C[-B^,^(0)]c} 

(C12) 

and 7(/3) is an error term arising from particles interacting accross the boundary of Bl^^^ (0). 
By the same argument as in 

flum, this term 

is negligible. Moreover, 

K ginp-m) 



m< E -^r( E e-* Hiv) 

m=L+l \e<S< m ) 

<[1 + (1)]^(5) E (P/») ro ( E e"** 00 Wwc**,,^)}) 



m=L+l r?g5( m ) 

where in the last inequality we use Lemmas IB.lHB.2l Now proceed as in (jB.19HTJ.20p , via the 
cluster expansion, to get 

K m 

<p((3)<l + o(l)]Cfi(S) EE E e-^ + ^ A -( A -^l 

m=L+Xj=X 2<fci kj<K (C.14) 

< [1 + o(l)] C e -0\PL+i-(ASp)] , 

where is the energy of a droplet with k particles that is closest to a square or quasi-square, 
r^+i = Hl+x + {L + 1)A, and the second inequality uses the isoperimetric inequality together 
with the fact that + kA is increasing in k for subcritical droplets. 
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On the other hand, by considering only those configurations in U > l that have a droplet 
with L + 1 paticles, we get 

<p(P) > [1 + o(l)] C /x(«S) e -/»tri+i-(A-^)] _ ( C15 ) 

Combining (jC. llj) and (|C.14lJC"T5j) . we conclude that linis-^ fXp{U > L) / [ip{S) = if and only 
if 

lim \Aa\ e-^^+i-^-^)) = 0. (C.16) 

B— >oo 



D Appendix: The critical droplet is the threshold 

In this appendix we show that our estimates on capacities imply that the average probability 
under the Gibbs measure fig of destroying a supercritical droplet and returning to a configu- 
ration in Sl is exponentially small in (5. We will give the proof for Kawasaki dynamics, the 
proof for Glauber dynamics being simpler. 

Pick M > £ c . Recall from (j!TTl) that ev M ,S L {o) = c/3(a)¥ a (t Sl < t Vm ) for a £ V M . Hence 
summing over a G &Dm, the internal boundary of Vm, we get using (|2.8[) that 

T,aedv M M^MgjPq (t 5l < rp M ) = CAP(S L ,V M ) ^ ^ 

Clearly, the left-hand side of (|D.1[) is the escape probability to Sl from <9£>ju averaged with 
respect to the canonical Gibbs measure pLp conditioned on dT>u weighted by the outgoing 
rate eg. To show that this quantity is small, it suffices to show that the denominator is large 
compared to the numerator. 

By Lemma [ 



CAP(S L ,V M ) < CAP(5 L , (S c \ C) UC + )=N \A g \ — e-^ (ig(S)[l + o(l)]. (D.2) 

On the other hand, note that 8T>m contains all configurations a for which there is an M x M 
droplet somewhere in A^, all L^-boxes not containing this droplet carry at most K particles, 
and there is a free particle somewhere in Ag. The last condition ensures that cg(a) > 1. 
Therefore we can use Lemma IB. II to estimate 

where Hj^2 is the energy of an M x M droplet. Combining (jTXMX3l) we find that the 
left-hand side of (|D.ip is bounded from above by 

A/3/ exp [-(3{H M 2 + AM 2 )\ 
which is exponentially small in j3 because T > H M 2 + AM 2 for all M > £ c . 
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